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Abstract 

In this paper we provide the theoretical basis for the problem of Griffiths- 
McCoy singularities close to the quantum critical point for magnetic ordering 
in U and Ce intermetallics. We show that the competition between Kondo 
effect and RKKY interaction can be expressed in Hamiltonian form and the 
dilution effect due to alloying leads to a quantum percolation problem driven 
by the number of magnetically compensated moments. We argue that the 
exhaustion paradox proposed by Nozieres is explained when the RKKY inter- 
action is taken into account. We revisited the one impurity and two impurity 
Kondo problem and showed that in the presence of particle-hole symmetry 
breaking operators the system flows to a line of fixed points characterized by 
coherent (cluster like) motion of the spins. Moreover, close to the quantum 
critical point clusters of magnetic atoms can quantum mechanically tunnel 
between different states either via the anisotropy of the RKKY interaction 
or by what we call the cluster Kondo effect. We calculate explicitly from 
the microscopic Hamiltonian the parameters which appear in all the response 
functions. We show that there is a maximum number N c of spins in the clus- 
ters such that above this number tunneling ceases to occur. These effects 
lead to a distribution of cluster Kondo temperatures which vanishes for finite 
clusters and therefore leads to strong magnetic response. Prom these results 
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we propose a dissipative quantum droplet model which describes the critical 
behavior of metallic magnetic systems. This model predicts that in the para- 
magnetic phase there is a crossover temperature T* above which Griffiths- 
McCoy like singularities with magnetic susceptibility, x(T) oc T~ 1+x , and 
specific heat Cy(T) oc T A with A < 1 appear. Below T*, however, a new 
regime dominated by dissipation occurs with divergences stronger than power 
law: X {T) oc l/(Tln(l/T)) and C V {T) oc l/ln 2 (l/T). We estimate that T* 
is exponentially small with N c . Our results should be applicable to a broad 
class of metallic magnetic systems which are described by the Kondo lattice 
Hamiltonian. 

PACS numbers:75.30.Mb, 74.80.-g, 71.55.Jv 
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I. INTRODUCTION 



In this paper we are interested in the quantum critical behavior of alloys of actinides 
and rare earths with metallic atoms. There is a large number of these alloys and they can 
be classified into two main groups accordingly to their chemical composition: 1) Kondo 
hole systems which are made out of the substitution of the rare earth or actinide (R) by 
a non-magnetic metallic atom (M) with a chemical formula Ri_a;M x (a typical example is 
Ui_ x Th 2: Pd2Al3); 2)Ligand systems where one of the metallic atoms (Ml) is replaced by 
another (M2) but the rare earths or actinides are not touched and thus have the formula 
R(Ml)i_ x .(M2) x (as, for instance, UCu5_ x Pd x .). In most of these alloys one usually has an 
ordered magnetic phase at x = which is destroyed at some critical value x = x* with a 
non-Fermi liquid (NFL) phase in the vicinity where strong deviations from the predictions 
of Landau's theory of the Fermi liquid are observed (see Fig.p. In this paper we propose 
that the origin for the NFL behavior is the existence of Griffiths-McCoy singularities at 
low temperatures. As we explain below, these singularities have their origin in a quantum 
percolation problem driven by the number of magnetically compensated moments. In this 
percolation problem magnetic clusters can tunnel in the presence of a metallic environment 
which produces dissipation. 

On the one hand, at Kondo hole system becomes an ordinary Fermi liquid with 

a temperature independent specific heat coefficient, j(T) = Cy/T, (by low temperature we 
mean T <C Ep where Ep is the Fermi energy of the metal), the magnetic susceptibility is 
paramagnetic and given by Pauli's expression, x(T) = xo, and the resistivity has the Fermi 
liquid form R(T) = Rq + AT 2 . On the other hand, at large x a ligand system is usually a 
heavy fermion [p]], that is, it can also be described by the Fermi liquid expressions but with 
coefficients 7 and Xo orders of magnitude larger than in ordinary metals. The nature of this 
heavy fermion behavior can be tracked down to the presence of the rare earths or actinides 
in the alloy. In the NFL phase the specific heat coefficient and the magnetic susceptibility 
show divergent behavior as the temperature is lowered. 
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The root for the understanding of the problem lies on the fact that in Landau liquids 
the thermodynamic and response functions are always well behaved functions of the tem- 
perature. This is clearly inconsistent with the behavior in the paramagnetic phase close 
to a quantum critical point (QCP). Exactly at the QCP divergences are expected since the 
system is ordering magnetically, thus, QCP can generate NFL behavior [^,[3] (albeit with fine 
tuning of the chemical composition). It turns out, however, that many times NFL behavior 
is observed away from the QCP and inside of the paramagnetic phase. One possibility is 
that the NFL behavior is due to single ion physics and therefore not at all related to the 
QCP physics. Another possibility is that there are still "traces" of the QCP inside of the 
paramagnetic phase. This is possible in the presence of disorder which can "pin" pieces of 
the ordered phase even in the absence of long range order. It is even conceivable that single 
ion physics and critical behavior co-exist close to a QCP. There is no clear consensus among 
researchers on the nature of this NFL behavior and the scientific debate is intense. The 
objective of this work is to shed light into this controversy. 

We organize the paper in the following fashion: Section II briefly reviews the problem of 
NFL behavior in U and Ce intermetallics; in Section III we write an effective Hamiltonian 
where Kondo effect and RKKY interaction appear explicitly and we discuss the problem 
of magnetic ordering and dilution within this Hamiltonian; Section IV contains a detailed 
discussion of the one impurity and two impurity Kondo problem and the extension of the 
problem to N spin clusters where we define a Kondo cluster effect; in Section V we propose 
the dissipative quantum droplet model which is the basis for the description of the problem 
of magnetic ordering in metallic magnetic alloys and we show that at low temperatures dis- 
sipation dominates the behavior of the system leading to universal logarithmic divergences, 
and that at higher temperatures non-universal power law behavior is expected; in Section 
VI we study the intermediate temperature range where quantum Griffiths singularities are 
expected; Section VII contains our conclusions. We also have included a few appendices 
with detailed calculations of the results used in the paper. 



4 



II. NFL BEHAVIOR IN U AND CE INTERMETALLICS 



The basis for the study of metals in the last 50 years has been Landau's theory of 
the interacting Fermi gas H. As consequences of Landau's theory the thermodynamic 
and response functions of the electron fluid are smooth functions of the temperature. The 
magnetic susceptibility has the form (we use units such that fi = ks = 1) 

X{ ' 4(1 + F °) ' 1 ' 

where g ~ 2 is the Lande factor for the electron, /ig is the Bohr magneton, Fq is an 
antisymmetric Landau parameter, 



m*kF 3p e 



is the density of states at the Fermi energy, m* is the effective mass of the quasiparticles 
(the effective mass is related to the bare mass, m, by a symmetric Landau parameter: 
m*/m= l + F-f/3), 

fc F =(37r 2 p e ) 1/3 , (2.3) 

is the Fermi momentum, p e is the number of electrons per unit of volume and Ep = 
hp/ (2m*). Moreover, the electronic specific heat, Cy, is given by the Fermi liquid expression 

^=^ = «- (2.4) 

Naturally these expressions resemble the ones obtained for the free electron gas where the 
bare mass m is replaced by the effective mass m* and Lande factor is renormalized by a 



factor of l/\/l + Fq. Furthermore, at low temperatures one expects the electronic resistivity 
to behave like 

R(T) =R + AT 2 (2.5) 

where Rq is the resistivity due to impurities and A is a constant. In a magnetic alloy where 
the magnetic moments are decoupled from the conduction electrons and do not interact 



directly among each other we expect an ordinary paramagnetic behavior for the susceptibility 
in the low field limit (T ^> gpsH where H is an applied magnetic field) 

X(T) = ^ (2.6) 

where C = PfS(S + l)(gyU B ) 2 /3 is the Curie constant (pj is the number of magnetic atoms 
per unit of volume, S is the magnetic angular momentum of the atom). At low enough 
temperatures (T -C gpsH) the susceptibility vanishes. 

The effect of disorder in ordinary Fermi liquids has been studied in great detail || 
and it was found that the results quoted above, especially the temperature dependences of 
the physical quantities, do not change appreciably (at least when disorder is weak enough 
to be treated in perturbation theory). Therefore, for a metallic system where Anderson 
localization J7[] does not play a role we still expect Fermi liquid behavior to be a robust 
feature. Thus, it is indeed very surprising that for such a broad class of U and Ce alloys 
(which clearly show three-dimensional behavior) deviations from Landau's theory are so 
abundant. 

In NFL systems it is usually observed that even in the paramagnetic phase the specific 
heat coefficient and the magnetic susceptibility do not saturate as expected from the Landau 
scenario. In all the cases studied so far the data for the susceptibility and specific heat have 
been fitted to weak power laws or logarithmic functions ||. The resistivity of the systems 
discussed here can be fitted with R(T) = Rq + AT a where a < 2. Neutron scattering 
experiments in UCus.^Pd^ || show that the imaginary part of the frequency dependent 
susceptibility, has power law behavior, that is, S(x(u>)) oc u 1 ^ x with A ~ 0.7, over 

a wide range of frequencies (for a Fermi liquid one expects A = 1). Moreover, consistent 
with this behavior the static magnetic susceptibility seems to diverge with T~ 1+A at low 
temperatures M. What is interesting about \JC\i5- x Pd x is that it has been shown in recent 



EXAFS experiments that this compound has a large amount of disorder |T0[ consistent with 
early NMR and /iSR experiments [TT| . 



From the magnetic point of view the U and Ce intermetallics show a rich variety of 
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magnetic ground states: antiferromagnetic as in the case of UCus-^Pda; for < x < 1 @, 
ferromagnetic as in the case of Ui_ 3; Th r Cu2Si2 for 0.6 < x < 0.85 |T^] and spin glass as 



in the case of \J\- X Y xP&z for 0.6 < x < 0.8 ||13|| . There are many indications that the 
magnetism in these alloys is inhomogeneous with strong coupling to the lattice. This seems 
to be the case of CeAl 3 f|TJ[] or CePd2Al 3 [I5| where /xSR experiments have shown clear signs 



of microscopic inhomogeneity. Moreover, because of the strong spin-orbit coupling in these 
systems a myriad of magnetic states with different ordering configurations are found [|l^ . 



On the theoretical side |16[ we can classify the approaches to NFL as of single ion 
character or cooperative behavior. Single impurity approaches for the NFL problem are very 
important because the Kondo effect [17j has been suggested as the source of heavy fermion 



behavior QXJ] in undiluted ligand systems (say, CeAl 3 ) [pl|. Nozieres and Blandin proposed the 



multichannel Kondo effect |19j as a possible source of NFL behavior. D. Cox proposed that 



such a multichannel effect of quadrupolar origin should exist in these systems |[20|| . Another 



source of NFL behavior based on single impurity physics is the so-called Kondo disorder 
approach which postulates an a priori broad distribution of single ion Kondo temperatures 
|TT| , PT|] . The cooperative approaches, on the other hand, stress the proximity to a critical 
point as a source of NFL behavior and was pioneered by Hertz [|2]|| and has been applied 
to many of the systems we discuss here [^-|26|. The problem of disorder close to the QCP 
has been discussed by many authors over the years. Hertz studied a classical XY model and 
showed that disorder is a relevant perturbation |27]] in dimensions smaller than 4. Bhatt 
and Fisher p8| and more recently Sachdev |29|] have argued that for Hubbard-like models 



an instability to a inhomogeneous phase should exist close to the QCP. A more recent work 
extending Hertz calculations to the case of disorder has reached similar conclusions, namely, 
that close to the QCP a new type of critical behavior, very different from the critical behavior 
of the clean system, should exist . 

It is clear from the experimental and theoretical point of view that magnetic interactions 
and the Kondo effect should be relevant for the physics of U and Ce intermetallics. More- 
over, it has been claimed for a long time that it is the competition between Kondo effect 



and magnetic order that determines the phase diagram of these systems. This idea was put 
forward by Doniach and we discuss it in more detail below. In trying to reconciliate 
the Doniach argument with the existence of disorder in these systems we proposed recently 
a new explanation for the NFL behavior in these systems which is based on the magnetic 
inhomogeneity due to this competition. In other words, if the Kondo effect is responsible 
for the destruction of long range order in the magnetically ordered system then at the QCP 
the competition between magnetic ordering and magnetic quenching is strongest [p2| . We 
have argued that close to the QCP, in the paramagnetic phase, finite clusters of magneti- 
cally ordered atoms can fluctuate quantum mechanically. In this case a so-called quantum 
Griffiths- McCoy singularities would be generated and zero temperature divergences of the 
physical quantities should be observed. In the next sections we explain exactly how this 
process occurs. 

Classical Griffiths singularities appear in the context of classical Ising models as due to 
the response of rare and large clusters to an external magnetic field [B3J . A classical model 



with columnar disorder was solved exactly by McCoy and Wu [33] in two dimensions and 
shows clearly the importance of Griffiths singularities close to the critical point for magnetic 
order. There are indeed experimental indications for the existence of such singularities in 



Feo.47Zno.53F2 and other Ising magnets ||35|| . Because of the relationship between the clas- 



sical statistical mechanics in d dimensions and the quantum statistical mechanics in d + 1 
dimensions the McCoy- Wu model can be mapped at T = into the transverse field Ising 
model which has been studied by D. Fisher |36|]. The transverse field Ising model is supposed 



to be applicable to insulating magnets since conduction electron degrees of freedom are not 
present. Senthil and Sachdev have studied the transverse field Ising model in higher dimen- 
sions on a percolating lattice and have found that quantum Griffiths singularities should be 
present in the paramagnetic phase in these systems close to the QCP ||37|| . Moreover, many 
physical quantities diverge at zero temperature as a result of the strong response of the large 
and rare clusters. Such anomalous behavior close to the QCP has been confirmed numeri- 



cally in dimensions higher than one |38[ and there is today strong experimental indication 



of such behavior in U and Ce intermetallics 
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III. EFFECTIVE HAMILTONIAN: RKKY x KONDO 

In our discussion of the physics of U and Ce intermetallics we are going to use the 
so-called Kondo lattice Hamiltonian which describes the exchange interactions between itin- 
erant electron spins and the localized magnetic moments: 

= £ eK (k)4 iK c k , K + £ £ J a>b (i)S a (i)cl(i)Tl K ,c K ,(i) (3.1) 

k,re i a,b,K,K.' 

where k — 1,2 labels the Spill st clt<3S, J a^b 

(i) are the effective exchange constants between 
the localized spins and conduction electrons at sites r,. S a (i) is the localized electron spin 
operator and J2 k ,k' c k(O t k,k' c k'(0 the conduction electron spin operator where c^ K annihi- 
lates an electron with momentum k and spin projection k (r b with a,b = x,y, z are Pauli 
spin matrices). Observe that we are allowing for the possibility of the conduction electron 
dispersion e K (k) to be dependent of the electron spin. Indeed, one can show that ( |3.1| ) can be 
obtained directly from the Anderson Hamiltonian when spin-orbit effects are included 
and the hybridization energy between localized and itinerant electrons is much smaller than 
the atomic energy scales [fil ]. One of the striking features of (|3.1|) is the fact that the ex- 



change interaction is not symmetric in spin space because of spin-orbit coupling. As we 
discuss below besides canted magnetism ( |3.1|) gives rise to what we call a canted Kondo 
effect. 

The main problem in studying the competition of Rudermann-Kittel-Kasuya-Yosida in- 
teraction (RKKY) and Kondo effect in the Hamiltonian Q3.1|) is related to the fact that 
both the RKKY and the Kondo effect have origin on the same magnetic coupling between 
spins and electrons. What allows us to treat this problem is the fact that the RKKY inter- 
action is perturbative in J/Ep while the Kondo effect is non-perturbative in this parameter 
and requires different techniques. Moreover, the RKKY interaction depends on electronic 
states deep inside of the Fermi sea, in addition to those at E F , while the Kondo effect is 
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purely a Fermi surface effect. Thus, it seems to be possible to use perturbation theory to 
treat the RKKY interaction. 

We will consider, for simplicity the case where the Fermi surfaces for the two species of 
electrons are spherical. The local electron operator can be written in momentum space as 

1 



CAl 



N k 



E e ' k ' rc k,- (3-2) 



We now separate the states in momentum space into three different regions of energy as 
shown in Fig.||], namely, Qq: kp ja — A < k < kp jCt + A; Q\: k < kp^ a — A; and f^: k > kp^ + A 
where A is an arbitrary cut-off. As we discuss below the value of the cut-off is related to the 
number of compensated spins in the system. But for the moment being we consider it as 
some arbitrary quantity we can vary, like in a renormalization group calculation. Observe 
that the sum in ( |3.2| ) can also be split into these three different regions. The problem we 
want to address is how the states in region Q Q close to the Fermi surface renormalize as one 
traces out high energy degrees of freedom which are present in regions Qi and f2 2 . We are 
going to do this calculation perturbatively in J/Ep ||43|| . For this purpose it is convenient to 
use a path integral representation for the problem and write the quantum partition function 
as 

Z = J DS(n,t)D^(r,t)Dij(r,t)exp{iS[S,^,tp]} (3.3) 

in terms of Grassman variables ip and i[) and where the path integral over the localized spins 
also contains the constraint that S 2 (n, t) = S(S + 1). The quantum action in ( |3.3| ) can be 
separated into three different pieces, S = S [S] + S [ip, ip) + «S/[S, ip, ip), where S = S [S] is 
the free actions of the spins (which can be written, for instance, in terms of spin coherent 
states 0), 

Sc#,^] = £ / 7rH^( k ^) (^ + /i-ea(k))5 Qi7 ^ 7 (k,^) (3.4) 
is the free action for the conduction electrons and 

<Sj[S, ip,ip} = J dtJ^Yl Ja,b{ri)S a (n, t)T b ai i) a (n, t)ij; 7 (n, t) (3.5) 



a, 7 n 
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is the exchange interaction between conduction electrons and localized moments. We now 
split the Grassman fields into the momentum shells defined above, that is, we rewrite the 
path integral as 

r 2 — — 

Z = / DS(n, t) J] D^(r, t)D^(r, t) exp {zS[S, ^}] } (3.6) 

where the indices 0,1,2 refer to the degrees of freedom which reside in the momentum 
regions Qq, Qi and Q2, respectively. The action of the problem can be rewritten as S = 
<So[S] + YU=o Soffit + "0) "0]- Notice the free part of the electron action is just a 
sum of three terms (essentially by definition since the non-interacting problem is diagonal in 
momentum space). Moreover, the exchange part mixes electrons in all three regions defined 
above: 

f 2 - 

Si = J2 dt ^a,fe(r ri )5 , a (r n ,t)r^ 7 ^ aii (r n ,t)?/' 7 ,j'(r n ,t) . (3.7) 

n J i,i'=0 

Since we are interested only on the physics close to the Fermi surface we trace out the 
fast electronic modes in the regions Qi and ^2 assuming that J a ^ fi. As we show in 
Appendix (0), besides the renormalization of the parameters in free action of the electrons 
in the region Qq, we get the RKKY interaction between localized moments. The effective 
action of the problem becomes: 

n a, 7, a, b 

+ J dtT^ b (r n -r m ,A)S a (r n ,t)S b (r m ,t) (3.8) 

n,m,a,b 

where T^ b (r n — r m , A) is the cut-off dependent RKKY interaction between the local moments 
and Jab( n ) i s the Kondo electron-spin coupling renormalized by the high energy degrees of 
freedom. As we show in Appendix this renormalization can be calculated exactly. For 
a spherical Fermi surface the exchange interaction between spins can be written as: 

r*,(r, A) = --— J2 Ja,cJb,cF(2k F r, A/k F )e~ r / e (3.9) 
hj F c 

where c is the number of magnetic moments per atom, n is the number of electrons per 
atom, I is the electron mean- free path and !F(x,y) is given in ( |A3 1|) . In Fig.|3] we compare 
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the usual RKKY with A = and the RKKY with finite A. At short distances, r <C kp , A -1 , 
the interaction is ferromagnetic and is renormalized by a factor of (1 — A/kp) 3 . Moreover, 
as shown in Fig.f| the first zero of the RKKY interaction is shifted from kpr m 2.2467 at 
A = to smaller values. At intermediate distances, that is, A -1 ^> r > kp 1 , the RKKY 
interaction decays like 1/r 3 as in the case of A = but the most striking result is that for 
large distances, r >> kp 1 , A -1 the RKKY interactions decays like 1/r 4 instead of the usual 
1/r 3 . Thus, for finite A the RKKY interaction has a shorter range. Another interesting 
result of a finite A is that at short and intermediate distances the RKKY oscillations are 
mostly antiferromagnetic. As one can see directly from Fig|| the ferromagnetic part of the 
RKKY is suppressed after the first zero. This could perhaps explain why most of the alloys 
of the type discussed here have antiferromagnetic ground states. 

We observe further that the perturbation theory here is well behaved and there are no 
infrared singularities in the perturbative expansion. Thus, the limit of A — > is well-defined. 
In this limit r^ 6 (r n — r m , A — > 0) becomes the usual RKKY interaction one would calculate 
by tracing all the energy shells of the problem. Observe that there are no retardation effects 
in tracing this high energy degrees of freedom since they are much faster than the electrons 
close to the Fermi surface and therefore adapt adiabatically to their motion. Hamiltonian 
( p.8[ ) is the basic starting point to our approach and contains the basic elements for the 
discussion of magnetic order in the system. Observe that the RKKY interaction depends on 
the electronic states far away from the Fermi surface while the Kondo interaction is a pure 
Fermi surface effect. While the RKKY interaction leads to order of the magnetic moments 
the Kondo coupling induces magnetic quenching. It is the interplay of these two interaction 
which leads to the physics we discuss here. The action (|3.8| ) has been used as a starting 
point for many theoretical discussions of rare earth alloys [jl"5| |. We stress, however, that the 
action ( |3.8|) describes only the low energy degrees of freedom of the problem and therefore 
the exchange constants that appear there can have strong renormalizations due to the high 
energy degrees of freedom. Moreover, as we explain below the cut-off A depends on the 
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number of compensated moments. 



A. Magnetically ordered phase: the role of RKKY 

The existence of local moments is not a sufficient condition for the existence of long 
range magnetic order. It is exactly the interaction between the spins which determines the 
ordering temperature T c of the material. There are various ways localized moments can 
interact: dipolar-dipolar interactions, kinetic exchange interactions due to the overlap of 
the f orbitals and RKKY interaction. Dipolar interactions are too small to account for the 
ordering temperature in these systems (they range from 100 K down to 10 K in the pure 
compounds) and the direct exchange between f orbitals is very weak since the spatial extent 
of the f orbitals is small (with the possible exception of the 5f orbitals of U). The RKKY 
interaction is by far the most important interactions in metallic rare earth alloys and it will 
be the only interaction we will consider in detail. 

As it is well-known and as shown in Appendix |^ in the ordered case (A = 0) the RKKY 
interaction decays like 1/r 3 and oscillates in real space with wave- vector 2/cp. The oscillatory 
terms have to do with the sharpness of the Fermi surface. Moreover, non-spherical Fermi 
surfaces will also lead to an angular dependence on the RKKY interaction with decaying 
rates which vary with the direction |j46| . The specific form of the RKKY interaction is 
not important in our discussion but the fact that the RKKY is an interaction which is 
perturbative in J//i and scales like N(0)J 2 . The exponential factor due to disorder was 
obtained originally by de Gennes |47| . 



Observe that the spin-orbit coupling generates an RKKY interaction which is anisotropic 
and therefore can give rise to canted magnetism. This effect is the analogue of the anisotropic 
spin exchange interaction, or Dzyaloshinsky-Moriya (DM) exchange interaction in in- 
sulating magnets which is obtained via the kinetic exchange between localized moments in 



the presence of spin-orbit coupling. The interaction ( |3.9| ) is an indirect exchange interaction 



in the presence of spin-orbit coupling. Thus, as in the case of DM interactions one expects 
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parasitic ferromagnetism within antiferromagnetic phases. This effect has been observed 



long ago in R-Cr03 systems p9| , p0| . The existence of ferromagnetic and antiferromagnetic 
coupling creates a very rich situation where many different magnetic phases are possible in 
the presence of a Lifshitz point ||51|| . Recent theoretical approaches for the NFL problem in 



CeCus-yAuy are based on the idea that the Lifshitz point in these systems is a QCP |23 
and therefore the quantum fluctuations associated with this point induce NFL behavior in 
the conduction band. 

The critical temperature of the system can be estimated directly from the mean field 



theory for (|3.1| ) and it is given by |52 



T c = 2S{S + 1} £ T M (R) cos (Q • R) (3.10) 

where Q is the ordering vector (Q = for ferromagnetism) and T^R) is the largest 
eigenvalue of r a ; ) (R). Observe that T c scales with T and therefore it is proportional to 
N(0)J 2 . This value of T c gives the order of magnitude of the transition temperature in 
Kondo hole or ligand systems, that is, the magnetically ordered Kondo lattice. In what 
follows we discuss the effect of disorder on the magnetic order in these systems. 

B. Magnetic dilution 

In Kondo hole systems the transition to the paramagnetic phase happens because the 
magnetic sublattice is diluted with non-magnetic atoms. In this case two main effects occur: 
(1) the magnetic system loses its magnetic atoms; (2) because the non-magnetic atoms do 
not have the same size of the magnetic ones there is a local lattice contraction or expansion. 
The first effect created by the dilution is to introduce disorder in the electronic environment 
and produce a finite scattering time r for the electron. As shown by de Gennes long ago 
||47|| , the RKKY interaction decays exponentially with the electron mean-free path, I = vft. 
Notice that this is only true if there is true magnetic long range order in the problem. In 
the case of a spin glass order this is argument is not correct . 
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The problem of destruction of magnetic order in a ligand system is more complex since 
the magnetic atoms are not replaced. In an alloy like UCu4_ a; Pd a ; the Cu atoms are replaced 
by the somewhat smaller Pd atoms. This difference between the Pd and the Cu leads to 
a local lattice contraction which modifies the local hybridization matrix elements. Since 
these matrix elements are exponentially sensitive to the overlap between different angular 
momentum orbitals one can have large local effects in the system. This change of local 
matrix elements induces changes in the exchange constants between the conduction band 
and the localized moments, J(rj), in ( p.l|) . 

Let us start with ( p.l|) and in the homogeneous ordered phase where we assume that 
Ja,b( T i) = J a ,b for all sites. The transition temperature is given by ( |3.10|) with T given by 



( p.9| ) with the electron mean free path determined by the extrinsic impurities in the system. 
As the Kondo lattice is doped, either by substitution of a magnetic atom by a non-magnetic 
one (as in the case of the Kondo hole systems) or just disordered (as in the case of ligand 
systems) the local coupling between the localized moments and electrons is changed from 
J a £ to a different value J a ^ (for simplicity we will assume just a binary distribution but this 
assumption can be easily generalized). In this case can rewrite ( |3.1| ) as H = H + H KL + H imp 
where: 

#o = X>K(k)ck, re Ck,K 
h kl = E E J a ,bS a (ri)cl(ri)T* jK ,c K ,(i) 

i a,b,K,K' 

Himp = ^Pi5J a)b S a (i)c ] K (i)T b KK ,c K i(i) (3.11) 

a,b,i 

where 5J a ,b = J a ,b — J a ,b where the summation over % includes all sites and p$ — 1 if a 
particular site i was changed by disorder and Pi = otherwise. 

Let us consider first the case where the system the ordered state is just slightly doped. 
In first order in x (the concentration) the problem reduces to a single ion problem. If 
\$J a ,b\ > \Ja,b\ then one has to treat first H + H imp and then add Hkl- It is obvious that we 
have a Kondo effect on the sites for which Pi = 1 with the original conduction band of the 
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system. In the case of lattice contraction we would have 5J a ,b > (J a ,b > J a ,b) and therefore 
an antiferromagnetic Kondo effect and a singlet is formed below a Kondo temperature Tk ~ 
Ep e~ 1 '( N ^' 5J > . All sites with Pi = I are magnetically quenched. If 5J a ,b < (J a ,b < J a ,b), 
which is the case of local lattice dilation, we would have a ferromagnetic Kondo effect and 
a local triplet state. Thus, there is an increase of the local magnetization of the system as a 
function of x. The effect here is very similar to the problem of enhancement of the magnetic 
moment by a highly polarizable metallic environment and creation of "giant moments" as 
it was discussed long ago by Jaccarino and Walker [54j and observed experimentally in 
Pdi^Nij and other similar alloys |55 . 



As it is well-known in a disordered system, the electron acquires a life-time, tk, due to 



the Kondo effect which, at zero temperature, reads ||17|| : 



J_ = 37t(SJ) 2 S(S+1)x 
t k 2E F 

where 6 J is the largest eigenvalue of 8J a ,b- This finite lifetime leads to a finite mean-free 
path, £k = vpTK for the conduction band motion. After the mean-free path is taken into 
account one proceeds as before (described in Appendix [AD to calculate the RKKY interaction 
which will be given by (|3.9|) with l/£ substituted by l/£ + 1/£k as in Matthiessen's rule. 

In the opposite limit of \5J a ,b\ < \J a ,b\ we have to consider H + Hkl first and then add 
to it H imp . This problem is more complicated because the ordered Kondo lattice problem 
has to be solved first. Here we just consider the simplest mean field theory in which the 
magnetic moments order along the z axis. The mean field Hamiltonian can be written as 
H = H M f + H imp where 

H M f = X! e «( k ) c k, K c k, K + l H s ( n iA ~ n i,l) + H s S z (ri)] (3.13) 

k,K i 

where H s = J z (S z {ri)) and H s = J z {{n^) — (n^)) are the molecular fields applied by the 
localized spins on the conduction electrons and by the conduction electrons on the localized 
spins, respectively. 

The problem described by Hamiltonian ( |3.13j ) can be easily solved for the case of ferro- 
magnetism as we show in Appendix (|B|). As a result the electronic degrees of freedom are 
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renormalized in different ways depending on geometry of the Fermi surface and the type 
of ordering one has. For ferromagnetism the main change in the problem is the change in 
the density of states for different electron species. In the case of antiferromagnetism the 
situation can be more complicated because one generates a well defined momentum Q and 
therefore Umklapp scattering is possible depending on the shape of the Fermi surface. A 
spin density wave state (SDW) can be generated and a gap can open on regions of the Fermi 
surface. Experimentally, the systems we are discussing here are metallic over the ordered 
phase which implies that the whole Fermi surface or perhaps large portions of it would re- 
main gapless. This is easily understood by the fact that the conditions for commensurability 
are hard to obtain in these systems which have rather complicated Fermi surfaces. There- 
fore, the conclusions we reach for the ferromagnetic case can be easily generalized to more 
complicated magnetic structures. 

Assuming that the system remains metallic in the magnetically ordered phase we see 
that Hk describes the Kondo effect on this new metallic band in the presence of a magnetic 
field. We show in Appendix [B] that at T = the local magnetic energies are H$ = —SJ Z and 
H s = SN(0)Jz/pf. Observe that while the magnetic field applied on the electron, H$, by the 
localized spin can be positive or negative depending if the local exchange is antiferromagnetic 
or ferromagnetic, the local field applied on the local spin, H s , by the conduction electrons is 
always ferromagnetic. In the paramagnetic case ((S z ) = ((n-j- — njj) = 0) the local state of 
the system is degenerate. That is, one has a quartet, made out of | | |,J|),| 1\1T) an d 

| t, JJ-) where f, j represents the conduction band spin and JJ-, the local moment spin. Since 
these are all eigenstates of S z and — their energies in term of the molecular fields are 



Ei 



SJ Z - 



SN(0)J* 



1A) 



E 2 



SJ Z + 



Pf 

SN(0)J* 



E 3 



SJ Z 



Pf 

SN(0)J% 



IA) 



Ei 



SJ Z + 



Pf 

SN(Q)J 2 Z 



(3.14) 



Pf 
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Observe that there is level crossing when J z = ±J C = ±p//iV(0). Furthermore, J c ~ cEp 
where c is the number of local moments per atom. Since J z <C Ep we can only have J z ~ J c 
when the density of local moments is very low (dilute limit). In general we expect J z <C J c 
in concentrated Kondo lattices. 

When J z > J c (antiferromagnetic coupling) the local state of the system is | |,Jj) which 
is separated from the state | L JJ-) by an energy amount 5E = 2SJ Z and therefore the Kondo 
effect is suppressed (notice that 5J < 5E and therefore the Kondo effect cannot bring 
these two states together). The problem is very similar to the usual ferromagnetic Kondo 
effect where quantum fluctuations are totally suppressed. Observe, however, that the local 
moment is compensated in the same way it would be if we just eliminate the local moment 
from the lattice. (Note that the electronic phase shift due to scattering by the impurity is 
zero in both cases.) Thus doping decreases the magnetization of the system as one would 
have in the usual dilution problem (the magnetic dilution of a ligand system is essentially 
identical to the dilution of the Kondo hole system). In the intermediate coupling regime 
of < J z < J c the two lowest energy states are | t?JJ-) an d | T,ff) which are separated in 
energy 5E = SN(0)J^/pj and no Kondo effect happens. The ferromagnetic case is somewhat 
similar with the difference that the local state of the system is | L Jj.) but also in this case the 
Kondo effect does not take place because the electron spin states are quenched by the local 
molecular fields. Therefore, in a ferromagnetically ordered lattice the Kondo effect is absent 
independent of the sign of the Kondo coupling. This state of affairs is very similar to the one 
found by Larkin and Mel'nikov in the case of magnetic impurities in nearly ferromagnetic 
Fermi liquids f56|| . In summary, we conclude that the case of \SJ a: b\ < \J a ,b\ there is no real 



Kondo effect in a magnetically ordered Kondo lattice. We note, however, that in the case 
of antiferromagnetic coupling the magnetization of the system drops because of a formation 
of | |, JJ-) states. 

We have seen that the ordered state of a Kondo lattice is destroyed via compensation 
of the magnetic moments either via the Kondo quenching or moment compensation. Thus, 
the percolation parameter in this problem is the density of quenched moments, pq. In 



percolation theory we assign a percolation parameter p which in the case of the Kondo lattice 
is essentially pq. Let us now consider the situation of one of the alloys mentioned previously 
where we chemically substitute the atoms by an amount x. In a Kondo hole system the 
number of magnetic moments decreases with the alloying because the magnetic moments are 
replaced by non-magnetic atoms. At the same time the number of compensated moments 
grows because of the changes in the local structure of the lattice and the increase in the 
Kondo coupling. At some particular value of x, pq reaches a maximum since it cannot grow 
beyond the actual number of magnetic moments which are left in the system. Moreover, at 
percolation threshold, p c , which is determined by the dilution and lattice changes, the last 
infinite cluster disappears and long range order is lost. Beyond this point only finite clusters 
can exist. Eventually both the number of magnetic atoms and the number of compensated 
atoms drop to zero. This situation is depicted in Fig.|5|(a). 

In a ligand system the density of magnetic moments is kept constant with chemical 
substitution because the magnetic atoms are not replaced. The number of compensated 
moments grows because of the local growth of the hybridization. At some critical value 
of pq we reach the percolation threshold p c and long range order is lost because the last 
infinite cluster of uncompensated moments disappears. Moreover, when pq grows beyond 
the threshold it will eventually reach the value pf and all the magnetic moments in the 
lattice are compensated. At this large value of doping one can find a heavy fermion ground 
state. We depict this situation in Fig.||(b). 

We have to be careful in interpreting the heavy fermion ground state in light of the 
Kondo model we are studying. As we mentioned previously, the dilution in a Kondo hole 
system leads to a trivial Fermi liquid state while in a ligand system it can eventually lead to 
a Heavy Fermion (HF) state which is not straightforward to describe. Since the HF state is 
non-magnetic it is clear that dilution can drive the system to such a state by increasing the 
local hybridization of the conduction band with the localized f-electrons. If the hybridization 
becomes of the order of the local atomic energy scales the Kondo Hamiltonian ( |3.1| ) is not a 
good starting point for the description of the magnetic correlations. We should work directly 
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with the Anderson Hamiltonian. This is usually the route taken by many approaches to the 



HF ground state [IS). The f-electrons mix with the conduction band electrons and the Fermi 
surface is large since it counts all the electrons in the system. In what we have discussed we 
considered only the Kondo Hamiltonian which does not contain this kind of physics since 
the occupation of the f-atomic states is fixed to be 1. Since we are interested mainly in the 
behavior of the system close to the magnetic ordered phase the Kondo Hamiltonian should 
give a good description of the problem but one has to be cautious about the transition from 
localized to itinerant behavior in these systems. 

Assuming that the Kondo lattice is a good starting problem we immediately see that the 
Hamiltonian generated by ( |3.8|) has the right properties associated with Doniach's famous 



argument |3l| , namely, there are two main energy scales described in ( |3.8|) : the RKKY 



strength T R and the Kondo temperature T# generated by J R . These two energy scales are 
local properties of the alloying procedure and they scale in a very different way with the 
bare exchange between conduction electrons and magnetic moments. The RKKY coupling 
is proportional to {N{Q)J{i)) 2 /Ep in the weak coupling regime while J R ~ J and therefore 
Tj((i) « Epe~ 1 /( N (°* )J ^' > (we are going to show below that this is slightly more com- 
plicated if we take anisotropy into account). If we plot these two quantities together as 
in Fig.|| we see that there is a critical value J c for the local exchange constant J(i) above 
which the Kondo temperature is larger than the RKKY energy scale and therefore as the 
the system is cooled from high temperatures the moment is locally compensated and the 
RKKY interaction for that particular place becomes irrelevant in the limit of large J{i)- For 
J(i) below J c the RKKY energy scale is larger and therefore as the system is cooled down 
the moment can locally order with its environment. (Note that for just two impurities, a 
partial Kondo effect occurs for finite J(i), and the ground state is always a singlet, even 
for ferromagnetic RKKY coupling p7| . Here, however, we are treating the case of a larger 
number of magnetic moments, tending toward a magnetic ground state for large RKKY.) 
We have to stress, however, that in the presence of disorder this effect is local and represents 
the a quantum percolation problem and has nothing to do with the homogeneous change 
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in the exchange J! Indeed, if one interprets the chemical alloying of the system as a simple 
change of the exchange over the entire lattice then the Doniach argument would predict 
an ordering temperature T/v as in Fig.^| which vanishes at a QCP where a transition from 
the ordered state to a fully Kondo compensated state happens. In this picture the QCP 
has nothing to do with a percolation problem where moments are compensated due to local 
effects. In our picture this is not possible since a local change in a coupling constant does 
not immediately imply a change of the "average coupling" constant. Thus, we believe that 
interpretations of the Doniach argument based on homogeneous changes in the couplings 
are actually erroneous. 

We also would like to point out that the same discussion carried out in terms of Hamil- 
tonian fl3.1|) can be done in terms of the effective action ( |3.8| ) by introducing a Hubbard- 
Stratonovich transformation and studying the saddle point equations by assuming a homo- 
geneous solution for the spin field. 

We have argued that by chemical substitution the order in a Kondo lattice can be taken 
away even when the substitution is not on the magnetic site because individual moments 
are compensated magnetically due to the distribution of exchange constants. At some finite 
concentration x* long range order is lost and the system enters a paramagnetic phase. Since 
the problem at hand is a percolative one, the paramagnetic phase can still contain clusters 
of atoms in a relatively ordered state. In the rest of the paper we are going to discuss exactly 
the physics of these clusters and how they respond to external probes. We are going to show 
how finite clusters give rise to Griffiths-McCoy singularities in these systems. 

C. The value of A 

In our previous discussion the cut-off A (cf. Fig. (2)) was considered as a free parameter 
but now we discuss this problem in more detail. It is clear from the above discussion that in 
the magnetically ordered phase we can set A = because all the electrons on the Fermi sea 
are participating in the order of the magnetic moments. As one starts to dilute the Kondo 
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lattice some particular sites with hybridization larger than average will be magnetically 
compensated by the Kondo effect. Therefore A has to increase in order to accommodate the 
electrons which participate in the Kondo coupling. Thus we expect A to increase with the 
number of compensated magnetic moments in the system. 

Let pq be the number of compensated moments per unit of volume at a given chemical 
concentration. As we have discussed previously this number is not given only by the Kondo 
effect alone but by the competition between the Kondo effect and the RKKY interaction in 
( p.8|) . If the electrons which participate in the Kondo effect are the ones in the region Oq in 
Fig.H then it is easy to see that 

d 3 k 



/d f£ 
——e(k F -A-k) 



(27T) 

(** - A) 3 m ,, 

= Pe-^^~ (3-15) 

which can inverted to give (using (|2.3| )) 

A=i_fi_^y /3 . (3.16) 

k F \ Pe J 

This equation gives a simple relationship between the number of compensated moments and 
the cut-off to be used in (|3.8j). When pq <C p e we have A/hp ~ pQ/(3p e ) <C 1 and A is 
very small compared with the Fermi momentum. Naturally, since each magnetic atom in 
the unit cell gives at least one electron to the conduction band, we must have pf < p e . This 
is true even when the electronic density changes as a function of the chemical substitution 
or percolation parameter p as shown in Fig.[5|. Obviously, we have pj > Pq and therefore 
Pq < Pe from which follows that A < kp. The extreme case of pq ~ p e (A « kp) indicates 
that Kondo process involves all electrons in the Fermi sea. Notice that the perturbative 
expansion in terms of J in the effective Hamiltonian ( |3.8| ) is still valid but the form of the 
RKKY interaction will not be the usual one, that is, ( |3.9| ), since it will depend strongly on 
A/ k f (as shown in Appendix |A]). 

Our discussion should be contrasted with the well-known exhaustion paradox proposed 



by Nozieres |58| : because the characteristic energy in the Kondo effect is the number of 
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electrons participating in the Kondo effect is supposed to be 

Pk„N{0) 3 Pe T K 

Pf Pf 2 Pf E F 

So, in order for all moments to be compensated one has to require this fraction to be 1. 
But because Tk <C Ef the condition in ( |3.17| ) is only observed at extremely small moment 
concentrations. As one increases pf there are not enough electrons at the Fermi surface to 
quench the magnetic moments and one should observe free magnetic moments. This paradox 
is based on an energetic argument which takes only the Kondo effect into account and is 
correct in the dilute limit. Indeed when pp > Pq — > we would have Tk ~ vfA. In the 
concentrated limit, where interactions start to play a role, it is misleading. The energy scales 
that produce A involve not only the Kondo coupling but also the RKKY interaction. Since 
the RKKY interaction involves states deep inside of the Fermi sea the energy scales involved 
in A can be rather large, of the order of the Fermi energy itself. Indeed, recent infrared 
conductivity measurements in YbInCu4 indicate that a large fraction of the Fermi sea must 
be involved in the Kondo quenching in clear contradiction to the exhaustion paradox and 
in agreement with our discussion [pH| . The calculation of pq thus involves a self-consistent 
calculation of the combined effect of Kondo and RKKY and goes beyond the scope of this 
paper. 



IV. 1,2,..., iV 

We have argued in the previous section that the local changes in the hybridization will 
lead to a finite density pg of compensated magnetic moments. The number of compen- 
sated moments depends strongly on lattice structure and the local changes in the electronic 
wavefunctions due to alloying. It is obvious, however, that alloying leads to a percolation 
problem and the number of magnetic moments drops as the ordered system moves towards 
the paramagnetic phase. Instead of working directly with pQ we define a percolation pa- 
rameter p. For p < p c magnetic order exists and for p > p c long range order is not possible. 
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p c is therefore the percolation threshold of the lattice. Thus, above p c there are only finite 
clusters of magnetic moments which are more coupled than the average. The probability of 
having N moments together is given in percolation theory by 



P(N,p)=N 1 -^f p (N/N i :) (4.1) 

where N% = {£,{p)/cl) d is the number of spins within a correlation length £, d is the spatial 
dimension of the lattice, is a percolation exponent (r 2 = 187/91 ~ 2.05, r 3 m 2.18), D is 
the fractal dimension of the cluster (jd = 1 + d/D). Close to percolation threshold p c the 
system is critical and therefore 

a \p-Pc\ 

where v is the correlation length exponent. The scaling function f p (x) is such that f p (x — > 
0) = 1 and for x ^> 1 one has 

f p (x) ss a;- ed ^ +Td e- c ^ Cp (4.3) 

where 2 , p >p c = 1, d 3 p>Pc = 3/2, 2tP<Pc = 5/4 and 6^ p<Pc = —1/9, c p is a constant of order 
of unit, Cp>p c — 1 an d Cp<p c — 1 — Observe that the exponential behavior given in ( fl.3p 
is the dominant part of the probability. 

The above equations are easy to understand in the dilute limit, that is, when p » p c . 
If c is the concentration of magnetic atoms in the systems then the probability of having iV 
atoms together is 

P(N,p » p c ) « c N = e ~ Nln{1/c) (4.4) 

and the probability is exponentially small. In particular, the probability of finding N nearest 
neighbor atoms in a lattice with coordination number Z is 

which reduces to ( [4.4]) when c — > 0. 
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It is also important to understand what happens close to percolation threshold (p — p c ~ 
+ ) where from ( |4. 1| ) and ( |4.3| ) we have 

f c N ) 

PWp^ocexpj-^^j. (4.6) 

Therefore, from Q4.4| ) and ( |4.(j| ), the exponential part of the probability can be written in a 
simple form 

P(N,p) (xe- K ? NCp (4.7) 

where 

k p » Pc « ln(l/c) 

>W « (C(p)/a)" D « b - Pel" - . (4.8) 

It is obvious from the above equations that the mean size of the clusters diverges at 
p = p c and that for p > p c the cluster size is determined by £(p). In this paper we are 
mainly interested in the physics of clusters in the paramagnetic case (p > p c ). For a given 
concentration c there are going to be clusters of all sizes but with different probabilities. For 
instance, for a cubic system [Z — 6), the number of isolated atoms (N = 0) becomes smaller 
than the number of dimers (N = 1) only when the concentration is larger than c > 0.14. 
In understanding the response of such an inhomogeneous to external probes we have to 
understand how a particular cluster responds to these probes. We are going to assume 
that the clusters only couple weakly to each other and that in first approximation they 
can be thought of isolated and permeated by a paramagnetic matrix. The weak coupling 
among the clusters, as we discuss at the end of the paper, can lead to glassy states which 
do not show strong deviations from a Fermi liquid ground state. In order to build up 
intuition about the physics of the clusters we discuss the case of 1 and 2 magnetic atoms in 
a paramagnetic matrix. After the discussion it will become quite clear how a cluster of N 
atoms should behave in the presence of a paramagnetic environment (this is the N impurity 
cluster Kondo effect). 
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A. 1 Magnetic Moment Kondo Effect 



The one impurity Kondo effect should occur in the case of Kondo hole systems at very 
low magnetic atom concentration and according to the argument given in Subsection |111 A 



also in the ligand systems in the ordered phase. Our discussion now will follow very closely 
the approach to the single impurity problem via bosonization techniques pT |. In order to do 



so we reduce ( |3.13| ) to an effective one dimensional problem via the bosonization technique. 

The bosonization procedure follows three different steps. In the first step we trace out 
the electrons far away from the Fermi surface exactly as in Section |T|. Since there is just 
one magnetic moment in the problem no RKKY term is generated and the only effect of this 
trace is to renormalize the Kondo couplings J. Indeed, as shown by Anderson et al. |)IJ for 
the one impurity Kondo problem the renormalization of the Kondo coupling is given by 

Jf = 8v F 5(J z ) (4.9) 

where 5{J Z ) is the phase shift of the electrons due to the scattering of a static impurity 
(corresponding to the Ising component of ( |3.8| )) which is given by 



:N(0)JA 



*5(J,) = arctan I wl - <■->-' \ (4 . 10) 



On a second stage we linearize the electron dispersion close to the Fermi surface: 

E Ka = E F + v F>a (\k\ - k F , a ) , (4.11) 

where we are considering the generic case where the spin flavors can have different Fermi 
velocities v FjCr . In order to linearize the problem we also have to introduce a momentum 
cut-off, A, which is a non-universal constant (anisotropies in the shape of the Fermi surface 
can be absorbed in A). The conduction band Hamiltonian is written as 

H c = J2 v F,aPcl a c p , a (4. 12) 

where c P)CT creates an electron with spin a, momentum |k| = p + k F and angular momentum 
I = 0. Moreover, since we are treating the problem of a single impurity it involves electrons 
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moving in a single direction associated with incoming or outgoing waves. Thus, in writing 
( |4.12| ) we have reduced the problem to an effective one-dimensional problem. In order to do 
it we introduce right, R, and left, L, moving electron operators 

1 f A 

^R(L)A X ) = -/jf= / dkc k±kFa , a e lkx (4.13) 

V Z7T J — A 



which are used to express the electron operator as 

Mx) = M^y kF - x + ^{xY~ lkF ' aX ■ ( 4 - 14 ) 

In any impurity problem the right and left moving operators produce a redundant description 
of the problem since they are actually equivalent to incoming or outgoing waves out of the 
impurity. Therefore we have two options: either we work with right and left movers in half of 
the line or we work in the full line but impose the condition ipR jCr {x) = ipL,a(~ x )- Following 
tradition we choose the latter approach. Thus, from now on we drop the symbol R from the 
problem and work with left movers only. The left mover fermion can be bosonized as 

Vv(x) = -Jw*-C*) (4.15) 
V ZTca 

where 

"M^) = J2\i^7 i( b p + ™ p )e™ px ~ h-c) (4.16) 



p>0 



pL 



and K a is a factor which preserves the correct commutation relations between electrons, 
that is, {ip a (x), ipa'(y)} = S(x —y)S<j,a> and the bosons obey canonical commutation relations 
[a k ,ap] = [6 fc ,6j] = 5 Pik . 

In what follows we are going to assume the simple case in which J a ^ = J a S a ,b where 
Jz < Jx < Jy In terms of the boson operators, the Kondo Hamiltonian ( |3.8| ) becomes 



H = v F ^p{a ] p a p + b ] p b p ) + 5v F ^p (b ] p a p + a+6 p ) + jfS z ^ \/-y( a fc + 

p>0 p>0 k>0 ' 71 L 



4) 



Ana 
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where Vp = ( v f,\+vf,i)/2 is the average Fermi velocity, 8vp = (v F ^ — v f,i)/2 is the mismatch 
between Fermi velocities in different spin branches, J + = (J x + J y )/2 and J_ = (J x — J y )/2. 
Notice that because we are allowing for different Fermi velocities the charge and spin degrees 
of freedom do not decouple from each other. Moreover, this Hamiltonian can be brought to 
a simpler form if one performs a unitary transformation 

(4.18) 



which transforms the Hamiltonian to 



H' = U l HU = v F J2 Pi a l a p + + Svf Eff (Mi a p + a lb. 



p>0 



p>0 



+ S 2 



( J? - nv F ) E ^(a* + 4) + Sv F E + b 



+ 



+ 



Ana 



fc>0 

J+ ( s + k\k } + s~k\k [ 

S~ K\K^e 4 ^v>°^ {ak ~ a * ) + h.c. 



4na 



(4.19) 



An important observation here is that the unitary transformation does not affect S z . The 
anti-commutation factors can be rewritten in terms of spin operators in which case we can 



rewrite [61| 



H' = v F ^2 P( a l a p + b l b p) + $v F ^2p {b ] p a p + a\b v 



p>0 



p>0 



2ira 



+ cr z 



jf - 7TV p 



E v -yfa* + 4) + Sv F E J -4(6* + b{) 



+ 



Aira 



(4.20) 



where a x = (ct + + a~)/2. Observe that ( [4 . 2 Q| ) describes the physics of a two level system 
coupled to a bosonic environment [ pT|j . 

When J z ^> 2v F and J_ = (the limit of large uniaxial anisotropy) we see from 
that Jjj, — > iiv F (Svp = 0) and the Hamiltonian reduces to 

J i 



H = 0~ x + V p E P a l a p 



(4.21) 



p>0 
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with the decoupling of the spin degrees of freedom to the bosonic modes (J x = J y = J + = 
J±). This is the dissipationless limit of the problem. Observe that in this limit the eigenstates 
of the system are eigenstates of a x , that is, the transverse field. 

The physical interpretation in terms of the Kondo problem is also very straightforward: 
in the limit of J z ^> 2vp the electron forms a virtual bound state with the localized spin. 
There are two states which are degenerate for an antiferromagnetic coupling 



1+) 
I") 



hi) 



(4.22) 



which are the eigenstates of o z . This degeneracy is lifted by the transverse field o x and one 
ends up with two states 



l«> = ^(l+>-|-» = ^(l^i>-|^t» 
W = ^(l+> + |->) = ^(l1U> + IM>) 



(4.23) 



which are the singlet and triplet states. 

Thus, we can define a tunneling splitting, A , between the two magnetic states and the 
coupling constant to the bosonic bath, a, which are given by 

J i 



a 



r T z ' 
1 Jr 


2 


, 2 




1 arctan 


7TVF- 




7T 



'nN(Q)J z 



(4.24) 



For simplicity we will consider the case of uniaxial symmetry in which J_ = and vp^ = Vp,[ 
so we rewrite ([4.20 ) as 



H' = A a x - nvpy/a a z \ -y( a k + 4) + v Fka\a k . 



(4.25) 



k>0 ' fc>0 

Observe that the magnetic moment flips at a rate given by I/Aq. The bosons with energy 
much larger than A (that is, with momentum close to A) will follow adiabatically the 
motion of the spin and their effect is to renormalize the fluctuation rate. The low energy 
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bosons (that is, with k pa 0) are too slow and cannot follow the motion of the spin. The spin 
can dissipate energy into this slow bosonic bath. In order to take into account the effect of 
the fast bosons into the motion of the spin we can consider the effects of bosons which live 
in a thin shell A — dA < k < A and treat their coupling to the spin in perturbation theory. 
It is simple to show that in this case the bare tunneling splitting A is renormalized to 



i 



A R = E c (^) " (4.26) 

where Ao = A(Ao) is the value of the tunneling at the bare value of the coupling constant 
and E c pa Ep is a high energy cut-off. The above discussion is valid at zero temperature. 
At finite temperatures the renormalization group flow has to stop at max(T, Ar) because 
after this point there are no bosonic modes to renormalize anymore. Thus, there is a 
crossover in the problem when the temperature becomes of order T pa Tr = Ar so that for 
temperatures larger than Tr the tunneling is suppressed and for temperatures smaller than 
Tr the tunneling is renormalized to Ar. This crossover temperature can be associated with 
the Kondo temperature as we discuss below. 

It turns out that the physics of ( [4.25|) is now understood and many physical quantities 



can be computed exactly For instance, it was found that, 



where 



{a z (t)) = (a z (0))e- lt cos(nt) (4.27) 



tv a ^ 2Tk ■ 2 ( na 
I a, An = sin 



Tr \2(l-a) 

p k . ( 
— sin 

Tr V I -a 



Q(a,A ) = — sin ( - ^ J 6(1/2 -a) (4.28) 



T R E c /A \Th; 



where 



K B T K = -Z = -Zl-^) (4.29) 

is the actual Kondo temperature of the system which depends on a non-universal cut-off 
energy scale E c . Observe that the ratio 
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only depends on a and it is universal. Moreover, it was shown that the zero temperature 
contribution of the magnetic moment to the susceptibility is simply p| 



Ximp (T = 0) = — L- = — L- (4.31) 
27rA^ 2itoi±k 



while the contribution to the specific heat is 

C imp (T) _ vr_a_ 7T 1 
T ~ 3 ~ 37V 

The imaginary part of the frequency dependent susceptibility is given by [65] 



(4.32) 



or ( w mY ■ I na \ u + (A QQ"\ 

y v imn u = sin — ^ tanhl — ) (4.33) 

LA pK n a \2{l-a)J (u 2 -Q 2 + T 2 ) 2 + 4T 2 Q 2 K 2 J K ' 

where 

(3 = 1/T. Notice that the Kramers- Kronig relation (x(0) = (1/ir) J du'V(l/uj')^s[x(w')]), 
leads to (|4.31|) . At high temperatures, T ^> T K , we have, for instance, 



Ximp(T) ~ 

C imp (T) oc ^ j . (4.34) 

The physics of the single impurity Kondo problem follows immediately from the exact 
solution ( |4.27|) : the electron spin acts as a transverse field on the local moment (see ( 4.25Q ) 



which flips (precesses) at a rate given by Q. In the process of flipping the magnetic moment 
produces particle-hole excitations in the Fermi sea which lead to an effective damping of 
the flipping process which is given by T. Observe that the flipping process is underdamped 
when a < 1/3 (since fl(a = 1/3) = T(a = 1/3)) and it is overdamped when 1/2 < a < 1/3. 
Oscillations cease to exist when a > 1/2. Going back to the Kondo problem we see that 
the oscillations can be classified as: underdamped if N(0)J z > 0.996, overdamped if 0.63 > 
N(0)J z > 0.996 and damped if N(0)J z < 0.63. Finally, just to make connection with the 
SU(2) Kondo problem let us observe that for J z , J± E c the Kondo temperature looks very 
similar to the SU(2) expression T K w £ c exp{-l/(iV(0)J)}. Indeed, from flOp we have 
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which is the form of the Kondo temperature for the anisotropic Kondo problem. Observe 
that the Kondo temperature of an anisotropic Kondo problem is not a single parameter 
quantity since it depends on the Ising component J z and the XY component given by J±. 



B. 2 Magnetic Moments Kondo Effect 

The single impurity Kondo problem discussed in the last subsection gives us a hint on how 
a cluster of coupled moments should behave, that is, one would expect renormalizations of 
the tunneling energy due to dressing of high energy particle-hole excitations and dissipation 
due to the low energy particle-hole excitations. It is however too simple because it does not 
contain the RKKY coupling between spins. The next level of complexity is the two impurity 
Kondo model |)7|] which has the first features of a cluster problem. 

We consider the problem of two magnetic atoms at a distance R from each other inter- 
acting via an RKKY interaction = T a 5 a ^ in the presence of a metallic host, that is, the 
so-called two impurity Kondo problem which is described by ( j3.8p : 

H d = J2^cl a c Ka + y £r*(R)S a (R/2)S a (-R/2) 

+ E ^[e^-* W ck^ (4-36) 

k,k',a 

where the couplings are renormalized by the trace over high energy degrees of freedom. 

In order to study this problem we reduce it to an one- dimensional problem by rewriting 
the electron operators as []S71 



1 /k-R\ iN , 1 . /k R N 

■ cos — - — + i{— ly sm — - — 



N e (k) \ 2 v ' N (k) 



c k , a (4.37) 



where j = 1,2 and Q is the solid angle and 



%W= li^f (4.38) 



32 



and {tp] ta {k), ipi,p(k')} = 2n5(k — k')5 a ,p8j,i- Furthermore, observe that all the momenta here 
are defined in a thin shell around the Fermi surface. Thus, we can linearize the band by 
writing = vp(kp — k) and rewrite the whole Hamiltonian close to the Fermi surface as 
(see Appendix [D]): 

d_ 

;.<> " ® X " J.I.O 



H = - iVF £ / dx^l a (x)—^, a (x) + £ rf s 1>a s 2>a + vJ2 4«(°KiMV 



2 



a,j,la,(3 



+ J-, a ^)rlK^iA^)Sa] (4.39) 

where r are Pauli matrices which act in the states j = 1,2, J s a are exchange constants, 
ipj tOC (x) is the Fourier transform of ( |4 . 3 7| ) and 



S a — $1,0, + 5*2^ 

5S a = S ha -S 2>a . (4.40) 

Moreover, we have explicitly introduced a new coupling constant V associated with impurity 
scattering and which breaks the particle-hole symmetry. This kind of term is unavoidable 
in any realistic model of impurities and plays an important role in what follows. 

Exactly as in the case of the single impurity given in ( |4.15| ) we bosonize the problem but 



take into account that besides f and j spin states we also have j = 1,2. This leads to four 
types of bosonic fields, namely pB[ ], 



<f) c (x) = (0i, T (a;) + <f>i,i{x) + </> 2 ,t(z) + 02,|(X>)/ 2 
<p a (x) = (0i,t(a?) - <Pi,l( x ) + 02, T (z) - h,i( x ))/ 2 
(f) f (x) = (0i, T (x) + (f)i,i(x) - faj(x) - h,i( x ))/ 2 

<Ps f (x) = (0 liT (x) - cj> u (x) - 2>T (x) + <f>2,i(x))/2 (4.41) 

which are associated with the charge, spin, flavor and spin-flavor currents of the problem. 
And again, like in the one impurity Kondo problem, we perform a rotation in the spin 
space in order to eliminate the spin currents. This can be accomplished with the unitary 
transformation 
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JJ = e -*(Si,*H-S2,«)*»(o) (4.42) 
in which case the Hamiltonian becomes 

f)th* \ 2 1 



h= v j_ y: jdx 

i=s,f,sf 



rr2/ \ , '-"Ti 

n, x + — 



2 

dx 



2V 

+ T ± (S ltX S 2!X + S 2>y S l!y ) + — cos($ s/ (0)) cos($/(0) - 9) 

v F J +:Z d§ s (0) v F J +>x 
H 5 ^ 2 H — cos($ fl/ (0))5 x 



+ — £-6S.-— — sm($ s/ (0))^ 



+ ^zii sin($ / (0) - 0) sin($ s/ (0))S 2 + cos($ / (0) - 9)S X (4.43) 

7ra 7ra 



where 



$i(x) = Ui(x) - [ dy Ui(y) 

\ J —oo 



9 = I CI dx ^ ($c(x) " $s/(x)) (4 - 44) 



and 



f r« 2vf (t „\ 
7ra 

J + ,z = J +iZ - 2tt (4.45) 

are the new renormalized couplings. Observe that ( [143] ) has strong similarities with its 
one impurity counterpart (|4.25|) but the bosons have not decoupled from the transverse 
field terms because of the appearance of bosonic modes $ s /. The most studied case of this 
Hamiltonian is associated with the NFL fixed point which is not our main interest since it 
does not describe the cluster physics we are looking for. The main problem with the NFL 
fixed point is that it is unstable to the particle-hole symmetry breaking operator V defined 
in ( |4.43| ). This operator is relevant under the renormalization group and drives the system 



away from the NFL fixed point. In the presence of scattering the system flows to a line 
of fixed points which represent the different phase shifts. This line of fixed points can be 
reached in different ways as we discuss below 
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Let us consider the case of V 7^ 0. As shown in refs. |66| V grows under the RG and 
accordingly to ( |4.43| ) the value of the <3> s / and $ / fields freeze at values such that 



$ s /(0) = 77777 

^(0) = im + e 



(4.46) 



where m and n are integers. By its definition in (|4.44|) we see that 9 is purely topological 
and depends only on the value of the fields at infinity, that is, depends on the boundary 
conditions. Thus it is clear that 9 is related to the phase shift the electrons acquire by 
scattering from the impurities. Therefore to each value of 9 we have a different fixed point 
and (f4.46 ) shows that in the absence of particle-hole symmetry the system flows to a line of 
fixed points as 9 varies in between and ir/2. Using (f4.46|) we find that Q4.43Q simplifies to: 



v F J x 



H — T z Si z S2, z + r_i_ (Si zS^ z + Si.yS^.jJ H 

7ra 



S r 



i=s,sf 



dx 



V F J +iZ d$ 8 (0) VFJm,zd$ sf (0) XQ 
H ^ k &z H ~ oo z 



(4.47) 



2n dx 2n dx 

where J± = J +: ± + J_ j_. Notice that the flavor degree of freedom decouples in this limit. 
This Hamiltonian has a form which is very similar to the two level system ( |4.25| ) we have 
studied previously. The main difference is that the spin operators for each impurity appear 
in a linear combination and therefore are not naively related to the two level system problem. 

The only trivial limit of the the Hamiltonian ( |4.47| ) is when J ^> T in which case the 
Kondo coupling has to be treated first. The impurities decouple and we essentially have two 



independent Kondo effects |p7| . We are interested in the opposite limit since in the cluster 
the RKKY is stronger than the Kondo effect. 

In order to understand the physics of Hamiltonian ( |4.47|) one needs to consider the type 
of RKKY interaction we are dealing with. Let us first consider the case where J +)Z = 
Jm,z — 0. The problem is purely magnetic since the bosons decouple from the spins. We 
can diagonalize the magnetic Hamiltonian since the it describes the simple problem of two 
interacting spins in a transverse field. When J± = the system decouples into two triplets, 
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I ft, ft) and | JJ-, JJ-) with S z = +1 and —1, respectively, and with energy T z /4; a S z = 
triplet, (| Jj, ft) + | ft, J|))/\/2, with energy -Y z /A + r ± /4; and a singlet 

I+,1> = ^(I1U>-|«», (4-48) 

with energy — r 2 /4 — Tj_/4. On the one hand, the transverse field J± splits the degeneracy 
of the 1 5^ | = 1 triplets down to 

|+,2) = -^(|^)-|1U>) (4.49) 
with energy T 2 /4 and (not normalized) 

|-2) = -±=(| i^A) + 1 1 1r» +e(l - I ^,1r» (4.50) 



with energy y (1^/4 — r_|_/4) 2 + J]_ + T±/A where J_l = v F J±/(2ira). The coefficient e is a 
normalization coefficient which depends on the relation between energy scales. On the other 
hand, while the singlet |+, 1) is still an eigenstate of the problem with energy — T z /A — T±/A 



and we get a new state 

~ y/2 



l ) = "4(1 t ^)) + e(| + \ M)) (4-51) 



with energy — y (f */4 — rj_/4) 2 + ± + T±/A. The splitting of the energy levels is shown 
in the diagram Fig. [5]. 

Let us consider first the case of antiferromagnetic coupling where T z > and large 
(T z ^> J mt x). The lowest energy levels are |+, 1) and |— , 1) which are separated in energy 
by -T ± /2 + 2Jl/f z « J 2 jf z (since T± « Ji/£ F and E F > f 2 ). Moreover, e « J ± /f z < 1 
and therefore the main effect of the transverse field is to split the degeneracy of the | ft, JJ-) 
and | Jj, ft) states. At low temperatures we just have to keep these two low lying states (we 
always assume T -C T z ) and introduce back the couplings J +j2 = J m>z . We see from Q4.471) 
that $ s decouples from the spins since 

S z \a, -a) = (S 1}Z + S 2tZ )\a,-a) = (4.52) 
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while the coupling J m>z gets a renormalization of a factor of 2 since 
SS z \a, -a) = (S ltZ - S 2>z )\a, -a) = 2a\a, -a) 



(4.53) 



which just tells that the sub-Hilbert space is spanned by the states | ff, Jj) and | Jj, f|~). The 
effective Hamiltonian can be written as 

/ o , \ 2" 



Haf = -y / <! •'• 



where 



U 2 sf (x) + 



dx 



v F J m . z d<S> sf (0) J\ 
2n ox r. 



(4.54) 





1,+) = 


+ |l,+> 


&x 


1,-) = 


HI,-) 


°z 


1,+) = 






1,-) = 





(4.55) 



Thus ( [4.54 ) has essentially the same physics as the initial Hamiltonian and its physical 
meaning is obvious, namely, it describes a dissipative two level system problem ( (4.25| ) where 
the two levels represent the two states | ij, Jj) and | JJ-, ft) of the magnetic moments, in order 
words, it describes a Kondo effect of the conduction band with the antiferromagnetic cluster 
made out of two local moments. Observe that the heat bath for the antiferromagnetic cluster 
is made out of the spin-flavor bosons. From the above Hamiltonian we immediately conclude 
that the Kondo temperature of the antiferromagnetic cluster is given by 



in. 1. 1 J'm,z j 



v 2 T 2 



a(J m ,z) \4Tr 2 a 2 T z E c 



n 



(27T) 



(N(0)J t 



R\2 



( sin(k F R) 
\ k F R 



(4.56) 



where we used ( P3|) . Observe that we have extracted a factor of 4 in front of a just to stress 
that a is coming from the coupling of the staggered moment of two impurities and therefore 
a scales like 2 2 a since y/a — > 2y / a in contrast with the single impurity case. 
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In the ferromagnetic case the situation is reversed since the two lying states due to the 
RKKY and transverse field are |2, +) and |2, — ) which are split from each other in energy 



by y 0V4 - r ± /4) 2 + Jl + r ± /4 - Y-JA » J|/r z . Moreover, we have 

S^cr, cr) = (S ljZ + 5 , 2, Z )|(T, cr) = 2\a,a) 
SS z \a,a) = (S 1)Z - S 2>z )\(T,cr) = 



(4.57) 



and therefore from ( 4.47 ) the $ s j fields decouple in the ferromagnetic case and the effective 
Hamiltonian in the sub-Hilbert space spanned by the triplet states | ft, ft) and | JJ-, JJ-) reads, 

/ o i \ 2" 



2 



dx 



where 



n i(x) + 



dx 



72 



o ?W + , z 0$ a (O) 

27T ax T 2 



(4.58) 



a x \2,+) = +|2,+) 




a x |2,-) = -|2,-) 




a 2 |2,+) = +|2,-> 




a,|2,-) = +|2,+). 


(4.59) 



Notice again that ( |4.58| ) describes a Kondo effect between the two low-lying states of a 
ferromagnetic cluster where the heat bath is provided $ s . The Kondo temperature of the 
ferromagnetic cluster is given by 



at(J. 



4n 2 a 2 f„E r 



4 1- 



N(0)J : 



(4.60) 



(2tt) 2 "V" 2 

where we used ( |4.45| ) and ( |L)3| ) . Notice the close resemblance of ( i.60|) and the single impu- 
rity problem where a is given by ( [4.241) . The only difference, like in the antiferromagnetic 
case, is a number in front which is associated with the number of spins involved. Indeed, 
the ferromagnetic problem is more directly related to the usual Kondo effect than the an- 
tiferromagnetic case. A direct comparison between the antiferromagnetic (|4.56j ) and the 
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ferromagnetic case ( f4.60| ) reveals a basic difference between the two problems: while dissipa- 



tion scales like (N(0)J z ) 2 <C f in the antiferromagnetic case it scales like 1 — N(0)J z « 1 in 
the ferromagnetic case. Thus, for the same coupling constants the ferromagnetic case has a 
stronger coupling to the bosonic bath. The antiferromagnetic, on the other hand, is weakly 
coupled and therefore dissipation is weaker. The same type of effect occurs in the problem 
of tunneling of magnetic grains where dissipation is more important for ferromagnetic than 
for antiferromagnetic granular systems H67|-|69|. 

The conclusion of the two-impurity Kondo problem is that the stable fixed points of the 
system represent either a ferromagnetic or an antiferromagnetic cluster undergoing a Kondo 
effect with the conduction band. The effective spin in this case is associated with the low 
lying doublet generated by the RKKY interaction which is split by the XY component of 
the Kondo coupling. Because the splitting of the doublet requires the flip of two spins the 
transverse field scales like the J]_ in direct contrast with the single impurity case where it 
is a linear function of the XY coupling. Moreover, the coupling to the bath is also modified 
since the flip of two spins leads to the production of more particle-hole excitations close 
to the Fermi surface. Although there is no way to define long range order for the two 
impurity problem it is clear that when the RKKY is dominant it is the "order parameter" , 
magnetization or staggered magnetization depending on the sign of the RKKY interaction, 
which couples to the heat bath. This trend seems to be easily generalized to more than 
two impurities, that is, although NFL fixed points are probably possible due to specific 
symmetries of the problem, when symmetries are broken a simpler Kondo effect of many 
moments is possible where the moments flip coherently in the spin field generated by the 
Kondo coupling. 

C. N Magnetic Moments - XYZ Magnetism 

Consider now the problem of N magnetic moments forming a cluster close to the QCP. 
If TV is large we can envisage this cluster as a large magnetic grain. The ground state of 
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the grain is just the classical one, that is, it is the fully ordered state for a ferromagnet 
(| fr, ff, ff) or | JJ-, JJ-, J|)) or the Need state for an antiferromagnet (| •fr, JJ-, -ft, J| ...) or 
I 4) fN JJ-5 fT •••))• Because of time reversal symmetry the ground state of a magnetic system 
has to be at least double degenerate. Like in the examples of the single impurity or two 
impurity the cluster can fluctuate quantum mechanically between the two degenerate states 
in the absence of an applied magnetic field which breaks explicitly the symmetry and bias 
one of the configurations. 

The tunneling between degenerate states can have many origins. In the preceding sec- 
tions we have discussed the tunneling due to the XY coupling of the Kondo interaction 
to the magnetic moment. Another more "mundane" source of interaction is the magnetic 
anisotropy in XYZ magnets. This kind of anisotropy exists even in insulating magnets and 
has to do with the interplay between spin-orbit and crystal field effects. To understand the 
origin of this anisotropy we can just look at the RKKY interaction alone in (|3.8|) . Observe 
that the RKKY interaction commutes with where = J2iLi S« is the total spin of the 
cluster. Thus, the eigenstates of the problem can be classified accordingly to the eigenstates 
of S 2 T , that is, S T (S T + 1) where S T = 0, 1, 2, N/2 if N is even or S T = 1/2, 3/2, .., N/2 
if N is odd. Because the cluster is in the ordered state it will select St accordingly to 
the interactions (St = N/2 for a ferromagnet and St = or 1/2 for an antiferromagnet). 
If the cluster is large it will behave like a magnetic grain and the total spin operator St 
behaves like a classical variable. Let us consider the simplest case of a ferromagnetic cluster 
(the antiferromagnetic case is essentially analogous) with iV atoms. Since the atoms are all 
locked together we can describe their spins as classical variables: 

S z ,i = S T sm(9) cos(0) 
S x ,i = S T sin(9) sm((f)) 

S y i = S T cos(fl) (4.61) 

where 9 is the angle the spins make with the Y axis and <fi is the angle in the X — Z 
plane. The RKKY interaction between the spins is given by (|3~8p with Y a ^ = T z 5 a ,b with 
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|r z | > \T X \ > \T y \ are the principal magnetic axis of the crystal. The energy due to the 
RKKY interaction can be rewritten in terms of the angles as 



E{9, <P) = N {T y - T x ) cos 2 {9) + (r 2 - r x ) sin 2 (#) cos 2 (0) + T 2 



(4.62) 



where T a = J2i^a(i) is the average exchange within the cluster. The energy ( |4.62| ) describes 
a magnet with Z easy axis and a X — Z easy plane if > T y > T x > T z . The energy of the 
cluster is minimized at = 0, it and 9 = tt/2 (that is, all the spins point along the Z axis). 
It is usual to rewrite fl4.62| ) as 



E(e,, 



N 



K ± cos 2 (9) - Ku sin 2 (9) cos 2 (0) 



(4.63) 



where K± = T y — T x > and Ku = T x — T z > are the anisotropy energies of the cluster 



(we have subtracted an unimportant constant from ( |4.63D ). 



The dynamics of a cluster described by the energy (|4.63|) is given by the well-known 



Landau-Lifshitz equations [70 



dS 



T 



dt 



S T x 



which, in terms of the angle variables, are 



5E 
5St 



dE 



(4.64) 



dt 
d9 



S T sin 9 89 
1 dE 



(4.65) 



dt St sin 9 d(p 

Of course, such equations do not allow for any tunneling. To study tunneling for such a 
magnetic grain we have to allow for solutions which are not classical in nature. The simplest 
way to do it is by the path integration method in imaginary time. The generating functional 



for the magnetic grain can be written as [59 



dfi[(j)]dfi[9]e 



S E 



(4.66) 



where Se is the Euclidean action 



S E = j dr ji£ T [l - cos(#)]^ + E(9, 



(4.67) 
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The tunneling process is now described as an instanton solution of the equations of motion 
for Se which interpolate between the minima of the potential. It can be shown that the 
tunneling energy is given by 



A^ = oJq\ cos(7rSV)| exp j— 2SVln 



1+KJK X + JKJK 



(4.68) 



where 



Lv = 2jK ll (K ll + K ± ) 



(4.69) 



is the attempt frequency of the cluster. Since St ~ N/2 we see that the tunneling splitting 
A can be written as 



A = u Q e 



(4.70) 



where 



7 = In 



1 + K\\/K ± + JK\\/K 



2 Ur, 



(4.71) 



Notice that the result ( |4.70[ ) is what one would expect from a WKB calculation for the 
tunneling splitting of N atoms. The factor of cos(ir St) has to do with the Kramer's theo- 
rem: if St is an integer (even number of spins) the magnetization can tunnel and split the 
degeneracy of the ground state; if St is a half-integer (odd number of electrons) tunneling is 
not allowed and the degeneracy of the ground state remains. This term has its origin in the 
first term in the Euclidean action Q4.67 ) and it is topological in origin [71]. Furthermore, 
the splitting is exponentially small in the number of spins in the cluster since St oc N for 
a ferromagnet. Notice that tunneling is only possible if K± ^ (that is, ^ T y ) which 
requires a magnetic cluster with very low spin isotropy, that is, XYZ magnetism. In an 
isotropic or uniaxial magnet (Heisenberg or XXZ, respectively) tunneling is suppressed. As 
we show below, only the Kondo effect can lead to tunneling. We also observe that at finite 
temperatures the system can be thermally activated from one minimum to another. The 
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process has the usual Ahrenius factor exp{— /3NKu} related to the jump of the system over 
an energy barrier of height NKu. Comparing this exponent with the one found in ( |4.68| ) 
we see that there is a temperature Ta above which thermal activation dominates and below 
which quantum tunneling dominates. This temperature is approximately given by 

NKu 

(4.72) 



2S T \n yi + K ll /K ± + ^K ll /K 1 

All the arguments presented here can be easily generalized to the case of antiferromagnetic 
grain |69]. 



In a metallic substrate dissipation due to particle-hole excitations in the conduction band 
plays an important role. Dissipation can be introduced in the problem via ( |3.8| ) f68f . Since 
within the cluster the moments are locked together we can rewrite ( |3.8| ) in Hamiltonian form 

as 

H = A A a x + e & c k ia ,o c k,a,o + Ja,bM a a a T b an ^ F n cl + ^ a c ka> o (4.73) 

k,a a,b k,q 

where J a ,b is the average exchange in the cluster, M a is the order parameter (S a (r) = 
M a a a cos(Q • r) where Q is the ordering vector in the cluster) and 

N 

F q = ^cos(Q-r i )e^ (4.74) 

i=i 

is the form factor of the cluster. 

The Hamiltonian (|4.73|) describes the Kondo scattering of the electrons by the cluster 



in the presence of a transverse field. As shown in ref. |]68| the transverse field is a relevant 
perturbation in the problem and all the processes associated with the XY component of 
the Kondo effect are irrelevant. In other words, in the presence of tunneling due to the 
anisotropy in the RKKY interaction the cluster flipping due to the Kondo effect does not 
play an important role. We have seen, however, that the splitting due to the RKKY interac- 
tion vanishes for Heisenberg or XXZ magnets and therefore the Kondo scattering becomes 
relevant. Thus, in the presence of we just have to keep the component of M a in the 
ordering direction (say, z). Because the Kondo coupling does not play a role any longer it 
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is possible to apply perturbation theory to Hamiltonian ( f4.73| ). It is easy to see that (f4.73|) 



maps again into the dissipative two level system ( f4.25| ) with A replaced by A^ and the 



dissipative constant a is given by a Fermi surface average pBfl : 

a = 2{N{Q)l z fMl JyJ^- cos(Q • r) cos(Q ■ r') 

= 2(iV(0)7,) 2 M^(5 Q , + 1)R {/ dre^^Mj ( 4 - 75 ) 

where the integrals are performed in the volume V c of the cluster. Thus, we introduce a 
cut-off term e~ r ^ L (so that V c = 8ttL 3 ) and perform the integration exactly: 

5c 



sm(kp |r — r'|) 
kp\r — r'| 



a = (N(0)J z ) 2 M? 



J Q,o 



1 + (2k F L) 2 

~ntn-n\(n -4- T] arrfan \(D — Oh „\ Tl\~ 

(4.76) 



1 / , . arctan[(Q + 2kp)L] arctan[(Q — 2k F )L] 

+ a tru2^ \ arctan(gL) 



AL 3 k 2 F Q \ v ^ ' 2 2 

The physical situation here is quite interesting. Although the XY coupling of the Kondo 
effect does not play a role, the anisotropy of the RKKY interaction introduces a transverse 
field. Thus, we can map back the two level system problem to a pure Kondo effect of the 
cluster! The only difference is that now in this "fake" Kondo effect the transverse coupling 
(what we called J± in ( [4.24[ )) is related to A^. The only difference from the original Kondo 
effect is that the transverse component is proportional to the original exchange as J 2 instead 
of J as in the true Kondo effect. Thus, one has a cluster Kondo effect with the RKKY 
coupling with a characteristic Kondo temperature 



a \E F , 

Observe that for a ferromagnet (Q = 0) the magnetization is M z = N/2 and in the large 
cluster size limit, kpL ^> 1, we find p8l 



a ^%k^^iT (m)1 ^ (478) 

where we assumed the cluster to be homogeneous with p/ = N/ (8irL 3 ) and p e is the electronic 
density given in ( ^.3| ). Observe that dissipation grows like a power law of the number of 
spins while the splitting decreases exponentially with the number of spins. 
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In the antiferromagnetic case Q ^ 0, the staggered magnetization is M z = N/2, and from 
( |4.76| ) we find that for large clusters (QL, 2kpL ^> 1) we have 



a AF (2i,N(l))J,) 2 (J^LA L 1 



+ C(l/L 3 ) (4.79) 



2 v ^' Q 2 -(2k F ) 2 QL_ 

Thus, in the antiferromagnetic case a has a singularity at Q = 2kp due to the Fermi surface 
effect. If Q < 2kp the leading order term is written 

a AF « ^(iV(0)7,) 2 (-^j N (4.80) 

and dissipation grows linearly with the number of atoms. Observe that the dissipation is 
substantially smaller than in the ferromagnetic case. For Q > 2kp we have 

/ l/3\ 4 ^2/3 

^.^(.(OPJ^J __ (4 . 81) 

which grows much slower with N as the previous cases. From this simple calculations we 
can conclude that the dissipation is much weaker effect in antiferromagnetic clusters. 
We can summarize these results as 

a = y (4.82) 

where 

^ = 4/3 

/ \ 1/2 , 

JV e = VfoF(^) — -J— (4.83) 



f f ) {N(0)J,)W 
for a ferromagnetic cluster, 

<p=l 

N = 4k2pQ U 84) 

7Vc n 2 p f (N(0)J z ) 2 

for an antiferromagnetic cluster with Q < 2kp and 

= 2/3 

(1 - (2fc F /g) 2 ) 3 / 2 / Q V 
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for an antiferromagnetic cluster with Q > 2kp. Notice that N c gives the critical number of 
spins in a given cluster above which the Kondo effect ceases to occur because for N > N c 
we find a > 1 which corresponds to the ferromagnetic Kondo coupling and therefore no 
Kondo effect. As we discussed before this case is related to the cessation of tunneling and 
the freezing of the cluster motion. Thus, the value of N c determines the largest size of a 
cluster which can still tunnel in the presence of a metallic environment when the anisotropy 
generated by the RKKY interaction is the source of quantum tunneling. 

We can estimate the value of N c for typical values of the constants. We assume Ep ~ 
10 A K, J z pa 300i^ ( which corresponds to an ordering temperature of the order T z pa 
J 2 JE F pa 1QK). Moreover, for U+^Cu^Pd we have p } w 1.2 x lCT 2 !" 3 and p e pa 
8.7 x 10~ 2 A -3 (which corresponds to kp pa 1.4A -1 ). Moreover, from neutron scattering we 
have Q ~ 0.8A -1 0. In the ferromagnetic case we find N c ~ 2 x 10 3 atoms and for the 
antiferromagnetic case N c pa 6 x 10 4 atoms! Thus, in the case of RKKY anisotropy a large 
number of atoms can quantum tunnel at low temperatures. We are going to see that the 
Kondo effect imposes much stronger restrictions on these numbers. 

D. N Magnetic Moments - XXZ and Heisenberg magnets 

We describe in this section that when the system has XXZ or Heisenberg symmetry (T x = 
T y ) the tunneling due to RKKY is suppressed. In this case the only source of tunneling, 
as we have discussed for the two impurity problem, is the Kondo effect itself. Again the 
XY component of the Kondo effect acts as a local magnetic field that flips the magnetic 
moment. This Kondo effect, in perfect analogy with the two- impurity Kondo problem, is 
due to the dissipative dynamics of states, say, |+) = | |,ft; |,ft; •••) and |— ) = 11)*^; T)4j •••) 
in the case of ferromagnetic coupling. The existence of quasi-degenerate low lying states 
separated from higher energy states is guaranteed by the fact that the cluster is effectively 
within the ordered phase and therefore states of the spins can be related by a finite number 
of spin flips. Now we can just borrow the results from the previous sections for the response 
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functions. In particular, the cluster Kondo temperature is given by 

where A (N) is the bare splitting between the two low lying states of the cluster and a(N) is 
the dissipative coupling to the bath. Since A Q (N) is the splitting between two states of the 
cluster where all the spins are flipped it has to scale like (Jj_/Tf) which corresponds to the 
energy required to flip iV spins. Furthermore, for a cluster the dissipation comes from the 
fact that the "order parameter" flips and produces a wake of particle-hole excitations close 
to the Fermi surface. It is clear that in this case in ( |5.4|) "order parameter" is an extensive 
quantity which implies that a(N) oc N 2 in complete accordance with the discussion of the 
two-impurity Kondo problem. Thus, we conclude that for a cluster one has 

f 



Ao(J\0 « r? (fj? ) =rfe-^ ln ( r ^) (4.87) 



z 



which by direct comparison with ( [4.70 ) we find 7 « ln(r^/Jj_) and ujq ~ Tf. Notice 



from Q4.45p that the RKKY interaction is strongly renormalized and in general Tf ^> J±. 
Furthermore, for the antiferromagnetic case the coupling of the cluster to the bath can be 
written as 



J, \ 2 / iV x 2 



aW ={ w< ) (4. 



while for a ferromagnetic cluster we expect (see (|4.60|) ) 

aW «*»(i -■£)' = (£)' (4 . 89) 



where N c ps Ep/J z (« (1 — -^) _1 ) for an antiferromagnetic (ferromagnetic) cluster. Observe 
that in a typical antiferromagnetic system N c w 100 — 1000 spins while for a ferromagnetic 
cluster N c w 1. We can immediately conclude that in the case of ferromagnetic clusters, like 
in the case where tunneling is generated by anisotropy of the RKKY interaction, dissipation 
will be so important that the cluster will freeze at low temperatures. We expect that ordinary 
superparamagnetism will occur - a situation much closer to a classical spin glass. Moreover, 
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observe that the Kondo temperature drops exponentially with the number of spins in the 
cluster and therefore for practical purposes the Kondo temperature drops to vanishing small 
values before the cluster size reaches N c . Thus, the region where dissipation occurs is rather 
narrow in the case of antiferromagnetic clusters. The Kondo effect, therefore, imposes much 
stronger requirements in the cluster size than the RKKY interaction. This is due to the 
fact that the Kondo effect requires the collective flipping of the cluster magnetization from 
electron scattering while the RKKY tunneling is a product of lattice anisotropy. Which 
effect is the dominant one in real systems depends on the lattice structure and the spin- 
orbit coupling. Furthermore, without dissipation the tunneling splitting vanishes only when 
the number of spins in the cluster diverges. In the presence of dissipation it vanishes for a 
finite number of spins given by iV c . 



V. THE DISSIPATIVE QUANTUM DROPLET MODEL 

We have seen in the previous section that the problem of a set of magnetic impurities 
interacting through a conduction band has two main features: tunneling and dissipation. In 
a simple mean field like picture there are three different energy scales in this problem: the 
ordering temperature which scales like T c oc N(0)J^, the tunneling energy which is given by 
Q and the damping energy which is given by V both given in ( |4.28| ). 



In order to understand how these energy scales affect the Kondo lattice consider a ligand 
system in the compositionally ordered phase. In this case the exchange between atoms is 
small enough so that the Kondo effect does not take place. As the system is doped with 
a metallic atom with size smaller than the original one the lattice locally contracts. The 
local matrix elements have exponentially large values and thus also the exchange J z . Then 
a particular region of the system can have an exchange parameter much larger than the 



average exchange in the lattice. This is the situation described by Hamiltonian (|3.13| ), for 
instance. If the exchange is locally very large then in face of the discussion of the single 
impurity Kondo problem it will be given by a Hamiltonian like ( [4.25| ) where the tunnelling 
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splitting is replaced by fi. Thus in the limit of large magnetic anisotropy the effective 
magnetic Hamiltonian which is generated from ( |3.8| ) is given by 

n,m n 

+ 5>*(r n ) (\ a (v n )b a + K(r n )bi) + $>J>1> Q (5.1) 

n,a a 

where a labels the relevant quantum numbers for the particle-hole excitations and A(r n ) is 
the local coupling constant of the spin to the electronic bath. When A = the bath decouples 
from the spins and the problem is mapped into the transverse field Ising model. We should 
stress at this point that ( |5.1|) can be only formally demonstrated for weak magnetic moment 
concentration. For larger concentrations it becomes an ansatz which has the correct features 
of the problems we are discussing. 

It is clear from the above that ( |5.1|) will drive the system through a quantum phase 
transition as Q increases. The simplest way to understand this effect is to rewrite the problem 
in path integral form and use the Suzuki- Troter trick to map the (i-dimensional quantum 
problem (|5.1| ) into ad+l classical problem by breaking the imaginary time direction into 
N T = (3/e pieces. The effective Hamiltonian associated with ( |5.1|) , after tracing out the 
bosons degree of freedom, reads: 

i=\n,m i=l n i^j=l "> \ l 3 ) 

where J n = ln[tanh(ef2 n )] and a n is the dissipation constant of each spin (observe that in 
the limit of e — > we have J n —> oo, thus, the mapping into the classical problem is valid 
when we take (3 — > oo). The original problem can now be thought as an anisotropic d+1 
classical Ising model with long range interactions in the imaginary time direction and short 
range interactions in the space direction. In the classical problem the sign of r n is irrelevant 
and for small J n the system will order magnetically (ferro, antiferro or spin-glass, depending 
on T n ) in d dimensions. Because the interactions are long-range in the imaginary direction 
the system orders in this direction as well. As in the classical case one assumes that the 
lowest excitation energy above the classical configuration is a droplet involving iV spins 
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which can be reversed at some energy cost E(N) which scales with the size of the droplet as 
N e . Moreover, we will assume that the droplets are very diluted and do not interact with 
each other in the spatial domain. It is clear that the problem is equivalent to independent 
droplets which are correlated in the imaginary time direction. The effective action for this 
problem is obtained directly from (|5.2|) by coarse-graining the spatial coordinate: 



Hn = E mN)S(ej) + K (N) S (ej) S (e(j + 1))] + ^ E ^f? ( 5 - 3 ) 

n [j=i 8 &k ^-3) j 

where S = — 1 corresponds to the ground state and S — +1 corresponds to the droplet 
excitation present at some imaginary time. The problem of a droplet at low energies is totally 
equivalent to a two level system problem. This approach is an extension of the quantum 
droplet model proposed by Thill and Huse |72| in the context of insulating magnets. We 
should point out that the quantum droplet model gives essentially the same results obtained 
numerically |38| and analytically |36| for the random transverse field Ising model. 

Notice that the energy barriers involved in the imaginary time are given by Aq(N) which 
is the surface free energy in d + 1 dimensions. It is obvious from the above discussion that 
in real time the dynamics of the problem is described by the Hamiltonian 

# = E \E(N)a z (N) + A (N)a x (N) - rr^{N)a z (N) E J^f{bk + 4) + E *4M (5.4) 

N \ k>0 ' nLl k>0 J 

which is the analogue of the two level system Hamiltonian (|4.25|) . Observe that a z is therefore 



related to magnetic order parameter of the system and a x is associated with the coherent 
flip of N spins within the droplet. In a ferro or antiferromagnetic state E(N) is proportional 
to the mean magnetic field of the spins which keep the droplet in the ordered phase. In a 
spin glass this energy is distributed accordingly to some distribution P(E(N)), say, 

P(E) = -l=e-^ 2 (5.5) 



where r = J 1^ m is the average interaction within the droplet. Moreover, as it is well-known 
in the context of the two-level system problem the bias E(N) is a relevant perturbation and 
will lead to ordering of the moments. In the rest of the paper we will interested only in the 
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paramagnetic phase where E(N) = and the droplets are just the clusters discussed in this 
paper. 

The statistical problem is now very similar to the one discussed in Section |V[ We 
have a distribution of energy scales (or cluster Kondo temperatures) which have their origin 
in a distribution of cluster sizes. Thus, contrary to the well-known distribution of Kondo 
temperatures approach to NFL behavior fTl~l , |2lH the effect is not of single impurity nature 
and the distribution is fixed by percolation theory. Using (|4.26| ), (14.871 ) and (|4.88|) we see 
that the renormalized tunneling splitting is given by 

/ 7iV + ln(W7u) ) \ 
l_(jv/iV c )V J 



A R (N) = We 



(5.6) 



where W ~ Ep. Usually it is not possible to invert the N as a function of A^ for a generic 
value of (p. Thus, for practical purposes the averages over the distributions of clusters have 
to be done with flO]) and flOp : 

N l-9 e -N/Ne 



P(N) 



(5.7) 



r(2 - e)N 2 f e ' 

We will consider first the example of the cluster Kondo effect for (p — 2 which serves as 
a good illustration of the effect of dissipation in the cluster problem. When p> = 2, ( |5.6|) can 
be inverted in order to give the size of the cluster in terms of the splitting or fluctuation 
time r = 1 / A R as 



N 



(iV c7 )2 +4\n(W/A R ) Huo/A R ) - N cl 



(5.8) 



N c 2\n(W/A R ) 

Notice that N vanishes when A# > ujq which corresponds to the smallest cluster (that is, 
N — 1 single Kondo impurity). Therefore the probability of finding a cluster with splitting 
A R is given by 

' dN 



P(A R ) = P(N(A R )) 



dA 



R 



oc 



k p N c {\n(W/A R ) \n(W/io ) + [y/(N c ^ + A\n(W/A R ) \n(uj /A R ) - N cl ] } 



x exp 



e-^ P )A R \n 2 (W/A R )^(N cl )* + 4\n(W/A R ) \n(u /A R ) 



-K 



p 2\n(W/A R ) 



(iV c7 )2 + 41n(W/A fl ) Hu /A R ) - N cl 



(5.9) 
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Observe that when N c — > oo we have 

P(A) oc A^- 1 (5.10) 

which gives a power law distribution for the energy levels. If Kp/7 < 1 the probability is 
divergent when A — > 0. As we are going to see in the next section it is this power law that 
gives rise to quantum Griffiths singularities. However, for finite N c the power law is modified 
at very small splittings. Indeed, when A# — > we find: 

P , A m .. N cKp (bi(W/<j ) + 7iVo) 1 ,,,,, 

P(A«-0)« ^^^^^ A^WA^) (5 ' U) 

which, unlike ( [5.10|) , diverges logarithmically. Thus, there is a crossover in the problem from 
the power law given by ( |5 . 1 ) for A# > A* to ( |5.11 ) for A# < A* where: 



A* = Jm7 exp |-I v /ln 2 (H//^ ) + ( 7 iV c )2| . (5.12) 

Thus, when N c — > 00 we have A* — > and the power law behavior dominates the entire 
range of energy scales. 



A. Magnetic properties: T < T* 

Since we now have the distribution of energy scales it is possible to calculate the physical 
properties of the clusters. For that we need the response functions for individual clusters 



which are given in Section [TV]. For the Kondo effect there are no analytic expressions for the 



physical properties which have to be usually calculated numerically |54|| . Here we will study 



only the asymptotic behavior. The temperature crossovers, however, can only be obtained 
numerically and in some special points which are not of practical interest. 

The magnetic susceptibility can be obtained directly from ( |4.31 ) and ( f4.34|) and for 
T < Ar(N) is given by 

* (T ' A * (iV))S3 2OT) (5 - 13) 
while for T ^> A#(A) we have 
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x (T,A r (N))k-^. (5.14) 



The average susceptibility is given by for T < u 

X(T)= / dNP(N) X (T,A R (N)) 



o 



N*(T) p(N) f N c p(N) 

dN — ^t4tt+ / dN-^r- (5.15) 
o 2ttA r {N) Jn*(t) AT v ' 

where A R {N*{T)) = T and we have approximated the integral by its two asymptotic pieces. 



Since A R < u in the high temperature limit (T > uj ) we have x(T) « 1/(4T) and therefore 
one has the usual Curie behavior. 

In what follows we will always assume that T <C Uq in which case N*(T) ~ iV c since 
A#(iV c ) = 0. Using ( |5.6p we easily find 

N*(T) knU- -r^ + ^WM ) (5 16) 



In this limit the second integral in (|5.15 ) can be computed immediately 



»° P(N) ^ N*- d (lN c + HW/u; ))e- N ^ D 

n*(t) AT ~ AT(2-6)(£/a) 2 - 9 T(<yN c + ip\n(pW)) { ' 

and therefore for T <^ T* = We~ lNc ^ we see that this term diverges like 1/ \T\n{W/T). 
The first integral in (|5.15 ) is more complicated but can be written in terms of exponential 



integrals and we find that in the limit of T <C T* it diverges like 1/ (T In 2 (W/T) and therefore 
is less singular than the second term. Thus, we conclude that when T<T* the susceptibility 
diverges at zero temperature as 

*< T > K T^WjT) <5 ' 18) 
which has stronger divergence a the power law singularity and it is independent on the value 
of <f. 

In the opposite limit, that is T ^> T* we see that N*(T) is very small and given by 

N*(T) ^-ln(^o) (5.19) 

7 
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which is what we expect from the usual Griffiths-McCoy case. We conclude therefore that 
power law behavior has to be observed at temperatures larger than T*. In systems where 
RKKY is the source of tunneling as we saw in Subsection [IV C| the value of N c is very large 
and therefore T* is effectively zero. Thus power law behavior will be predominant over all 
realistic low temperatures. In systems where the Kondo effect is responsible for tunneling, as 



described in Subsection |IV D| we can have N c pa 10 if J z ~ 1, 000 K and 7^2 (for ordering 
temperatures of the order of pa 10 K and J± ~ 1 K) and therefore T* pa 0.5 K which is 
a reasonable energy scale. Notice, however, that T* is an exponential function of N c and 
therefore a change to J z ~ 500 K (which is probably a better estimative for the exchange in 
rare earths) leads to T* pa 10~ 9 K\ 

B. Specific Heat: T < T* 

Once again analytic expressions for the specific heat of a Kondo system are not known. 



We will use the asymptotic forms given in (|4.32|) and (4.34). For T <C A R (N) we have |64 



. C V (T,N) 7c a(N) 
7 T, N) = VV ' 7 pa - - \ \ 5.20 



while for T > A R (N) 

<y(T_ N) Pa , 

where we used Q4.26 ). Exactly as in the case of the susceptibility we break the integral of 



/A \ 2 / TP \ 2 ( 1 "«( Ar )) 

7(T,JV)» (5.21) 



the specific heat into two pieces: 

n rN* (t) P(N)a(N) /An\ 2 r N ° /b\ 2 MW) 
7(T) - I / dN^P^l + ^ ) dNP(N) ( ^ ) . (5.22) 



3 70 A R (N) \EJ Jn*(t) v y VT 

For T <ti T* the integrals can be evaluated as in the case of the susceptibility. The main 
difference here is that the second integral in ( |5.22| ) which gave the most divergent contribu- 
tion in the case of the susceptibility diverges only like l/ln(/3W / ) because a(iV c ) = 1. Thus 
in the case of the specific heat the dominating term is the first one which diverges like 
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which is a weaker divergence than the susceptibility. In the limit of T ^> T* we recover the 
power law behavior characteristic of Griffiths-McCoy singularities. Other physical quantities 
can be directly calculated from the probability distribution ( |5.7| ) and the dependence of 
Ar(N) on N given in 



VI. T > T*: QUANTUM GRIFFITHS SINGULARITIES 

We have shown in the previous section that at temperatures smaller than T* the cluster 
dynamics is essentially dissipative and the singularities in the thermodynamic functions have 
logarithmic character. A reasonable estimate of T* gives a very low temperature which in 
many systems is experimentally not reachable. Thus, for most systems the physics in a 
wide temperature range from T* to T « T where T is the average RKKY interaction in the 
system the physics is non-dissipative. For temperatures above T the cluster does not exist 
as a well-defined object since temperature fluctuations can excite isolated spins within the 
cluster. In the temperature range T* < T < T we can effectively treat the system as non- 
dissipative, that is, set a = (or equivalently, N c — > oo). In this crossover temperature scale 
the physics is highly non-universal and depends strongly from the distance from percolation 
threshold. Indeed, consider the probability of finding a cluster with an energy splitting A 
which is given in ( |5.1(J| ). Since there is a broad distribution of energy scales in disordered 



systems it is often convenient to use a logarithmic scale so that 

P(ln(A)) oc A Ap (6.1) 

where 

A P = ^ (6.2) 

7 

is an important exponent of the theory. Equation ( |6.1|) is quite revealing if one considers 
the probability of finding a spin in a cluster of size N or excitation energy A. Since the 
probability is proportional to the size of the cluster one sees that this probability is 
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P L (ln(A)) oc iVA K f /7 oc L D A K ? h 

oc (LA 1 ^) (6.3) 



where 



Z P = £ (6-4) 

is the so-called dynamical exponent. Observe that Z p depends directly on X p which is a 
non-universal quantity which depends on many microscopic details as one can easily see 
from the definition ( |6.2|) . Thus Z p varies over the phase diagram. Indeed, using ( |4.8|) and 

n 



7 ~ 

Z/ P>Pc ~ 



7D 



ln(l/c) 

Zp~p. « lD{ap)/a) D » b _^ D - 00 (6.5) 

which diverges at percolation threshold. The physical meaning of the dynamical exponent is 
clear from ( |6.3[ ): it defines the relationship between length and energy scales in a correlated 
volume L of the system. The thermodynamic and response functions of the system are 
directly related to Z p . Finally, we have to point out that the dynamical exponent in the 
Griffiths phase has a different meaning than the one usually used for clean systems @ (thus 
we use the capital letter Z p for the dynamical exponent in the Griffiths phase). We should 
stress once again that the behavior of the system for T < T* also corresponds to Z — > 00 
but unlike the non-dissipative system Z is divergent away from the critical point since it is 
independent of k p which vanishes at p — p c . Observe that even at this point ( |5.11|) is finite. 



All the thermodynamic functions can be now calculated from the knowledge of A p or 
Z p . Although we were able to find microscopic expressions for 7 and ujq the relationship 
between p and Z p depends on details of the percolation problem. Thus, another way to face 
the problem here is to assume X p as a phenomenological parameter which can be varied as 
doping is varied and which vanishes when p — > p c . Thus, from now on we drop the subscript 
p on A p . In the rest of the section we give an explicit calculation for the response functions 
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in the crossover temperature range. Using ( |5.10|) and ( |4.70|) we find that the normalized 
distribution function ( |5.10| ) is given by 



[ln^o/A)] 1 " / A\ A-1 

p(A) = r[2-^]( 7 iv^ U e ^" A )- 



(6.6) 



Notice that the Heavyside step function appears because the largest value for the tunneling 
splitting is ojq. 

The calculation of physical quantities is now rather simple because we have just a two 
level system Hamiltonian to deal with. For generality let us consider the problem of a 
magnetic field H applied along the easy axis of the cluster so that the effective cluster 
Hamiltonian can be written as 

H c = E H a z + Aa x (6.7) 

where Eh is the magnetic energy of the cluster. This magnetic energy can be written as 
Eh = m c H where m c is the mean magnetic moment of the clusters and it depends on the 
magnetic nature of the cluster and the number of spins in the cluster. For a ferromagnetic 
cluster we would have m c = qhbSt- If the cluster is fully aligned then St = N/2. For an 
antiferromagnetic cluster in a percolating lattice we can have the following possibilities: (1) 
the cluster has random fluctuations of the spins (uncompensated) in its surface in which 
case the net magnetization is proportional to (N) l l D (this result can be obtained from the 
fact that the average size R of cluster scales like R oc N 1 ^ and therefore its area scales like 
A oc N 2 I D leading to fluctuations or the order \[A oc N l / D as iV ^> 1); (2) the cluster has 
uncompensated spins in its volume so that the net magnetization proportional to (iV) 1 / 2 . 
In general we write m c = fiBqN^ where q gives the magnitude of the average moment in 
the cluster (for a spin glass it is the Edwards-Anderson parameter) and <fi = 1(1/2) for 
ferromagnets (antiferromagnets and spin glasses). Thus, using ( [4.70| ) we can write 



E H (A) = qnB 



1 



H. (6.8) 



7 V A. 

Notice, however, that the dependence of Eh on A is rather weak and since A is cut-off from 
above by c^o, H or T, we can safely replace the expression above by 
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E H (H,T)^qii B 



In 



H . 



7 \max(H } T) 

In doing so we disregard higher logarithmic corrections in the expressions given above. 
The cluster magnetization is obtained from ( |6.7| ): 



(6.9) 



M(A, H, T) 



El + A 2 



tanh (pyjEjj + A 2 ) 



(6.10) 



the static magnetic susceptibility can obtained directly from ( |6.10| ): 



Xo {A,H,T) = fi 3 Bl 



PE 



H 



sech z (3JE 2 H + A 2 + 



A 2 



tanh pJE 2 H + A 2 



E 2 H + A 2 V v H J (E 2 H + A 2 ) 3 / 2 
the imaginary part of the frequency dependent susceptibility for uo > is (see Appendix 
(0)): 



(6.11) 



S\x{u>,A,H,T)] = 2n- 



A 2 



(E 2 H + A 2 ) 
and the specific heat is given by: 



tanh (3JE 2 H + A 2 5(u - 2JE 2 H + A 



(6.12) 



CV(A, H, T) = f3\E 2 H + A 2 )sech 2 f3JE 2 H + A 



(6.13) 



Any physical constant can now be calculated directly from the expressions above and the 
distribution of tunneling splittings given by (IT 



1. Magnetic Response 



From ( |6.10|) and ( |6.6| ) we see that the average magnetization is given by 



r^o A A_1 

M(H,T) ocH dA 

'o JE 2 +A 2 



tanh ( ^E 2 H + A 2 ) [In ^ 



oc H^- x I 
J0 



dx- 



tanh(x) 



:[ln 



(3(jjQ 



n 1—6 



'Wh \x 2 -{(3E H ) 2 f- x/2V ^ x 2 - (pE H )< 

Observe that the magnetization has a scaling form 



M(H, T) 



H (H T 
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(6.14) 



(6.15) 



where f\(x,y) is a simple scaling function. 

Let us consider first the high temperature case (T ^> u , E H ) in which case the integral 
simplifies to 

M m (6.16) 

which is just the usual Curie behavior. Observe that the high temperature results have to 
be understood with a little grain of salt since the cluster decomposes when the temperature 
becomes of the order of the characteristic coupling between spins (which is of order of order- 
ing temperature of the pure system, T c ). At high temperatures one has indeed paramagnetic 
behavior like the one described by ( |6.16|) but the Curie constant can be quite different from 



the one we would get from the cluster calculation since it comes from individual atoms. 
Thus, high temperature behavior in the cluster picture only makes sense if T c ^> T ^> u 
which puts a strong constraint on the temperature range. Moreover, at high temperatures 
the clusters can be thermally activated over the barrier and the two level system description 
is not complete. We therefore focus on the case where uq 3> T, Eh- Notice that in this case 
the scaling function in fl6.15| ) depends only on the ratio of H/T. 

Let us consider the case where u ^> T, E H . If T E H (low field limit) ( |6.14[ ) can be 
safely replaced by 

Vx iM^ H (617) 

Observe that the magnetization is linear in the field and diverges at low temperatures as 
rpx-i f or ^ < x. Thus, the scaling function is f\(x, 0) ~ 1 when x — > 0. 

On the other hand, if ujq 3> Eh 3> T (high field limit) the integral has a different 
behavior, namely, 

1-0+4 



POO 

M(xH(P) x - x dx 

J0E„ 



\n^J^ 

y/x2-(J3E H )2 



2 _ frtE?„Wl-V 2 



x^fln^) (G.1S) 



PE H (x 2 - (J3E B 
1-0+4, 
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which shows that the magnetization scales like H x and therefore the susceptibility behaves 
like if A_1 . Thus, the scaling function behaves like f\(x,0) —>■ x A_1 when x — > oo. The 
crossover from ( |6.18| ) to ( |6 . 1 7| ) occurs when Eh ~ T. 



The imaginary part of the average frequency dependent susceptibility can be easily calcu- 
lated from fl6.6| ) and fl6.12| ) since the only contribution comes from the Dirac delta function: 

t nh (^ w \ ( \ 
Im[ X (u,E H ,T)} oc U - (2E H ) 2 ) X/2 ^LiZZ (in 2u ° =) Q(uj - 2E H ) . (6.19) 
v 1 u \ ^-(2E H y) 

This function can be measured in momentum integrated neutron scattering experiments 

since it is the local response function. Observe that for Eh = it reduces to 

tanhf^r) / 2un\ 1 - 9+ ^ 
Im[x(u,0,T)} oc i—x ' — ) (6-20) 



which is rather simple function. The fact that this function vanishes for u < 2Eh has to do 
with the blocking of the tunneling by a magnetic field which biases one specific magnetic 
configuration. Observe that the application of a magnetic field suppresses the quantum 
fluctuations and therefore destroys the NFL behavior caused by the clusters. 

Another important quantity for the characterization of a Griffiths singularity is the non- 
linear magnetic susceptibility which can be calculated from ( |6.15 ) as the third derivative of 



the magnetization. Due to the scaling behavior it is obvious that the leading behavior is 

XsCT) oc JL (6.21) 

which therefore diverges more strongly than the linear susceptibility. 

In certain experiments, like the Knight shift K(T) measurement in NMR, it is possible 
to observe the width of the distribution of energy scales by studying the distribution of local 
susceptibilities in the system. Indeed, we can write: 

«i> = o s Jsn (6 . 22) 

X K 

where D is a constant which depends on the range of the interactions. For long range inter- 



actions D = 1 while for short range interactions one has D > 1. $x(T) = yx 2 (T) — x(T) 2 
can be calculated directly from ( |6.11| ) and ( |6.6|) 
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for T — > 0. This result shows that as the temperature is lowered the width of the distribution 
of susceptibilities grows and eventually diverges at T = 0. In fact any moment of the 
susceptibility can be calculated from ( |6.11| ) and ( |6.6|) . 



2. Neutron Scattering 

Consider now the problem of scattering of neutrons by the magnetic clusters discussed 
here. The differential cross-section for neutron scattering is given by the usual expression 

n(^)'Ui-) (6.24) 



dVtdu \mNC 2 ) k 

where q = k — k' is the momentum transfer and u = {{k') 2 — k 2 ) / (2m n) the energy transfer 
in the scattering. S(q, u) is the dynamical form factor which can be written as 

S( q ,u) = E J dte^^<*^M(P)Sj(t)) , (6.25) 

where F(q) is the nuclear form factor. It is usual to rewrite the spin-spin correlation function 
as (Si(0)Sj(t)) = ((Si(O)Sj(t)) — (Si)(Sj)) + (Si)(Sj) so that we can divide the dynamical 
form factor into a static and a dynamic part: 5(q, uS) = Ss(q)<5(u;) + S^q, uS) where 

S s (q) = —pL'£ i e i *to-*A(S t )(S j ) (6.26) 

i,3 

gives the static response. We are only interested in the dynamic response which can be 
rewritten in terms of the imaginary part of the susceptibility as 



S D (q,u) 



2F 2 (q 

Ntt 



1 



3{ X (q,^)}. (6.27) 



1 _ e -Pu>_ 

In the previous subsections we have calculated the q = component of S{x(q, u))} but 
for neutron scattering we need a more detailed description. In order to do that we again 
consider the clusters as independent scatters. Moreover, since the cluster scatters as a whole 
we write for a cluster of iV atoms that 
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9 {;ov(q, u)} = F 2 N (q)% {X(w, A(iV), T)} (6.28) 
where F^(q) is the form factor of the cluster and S {xi^i A(iV), T)} is given in Q6.12 ). The 



main theoretical problem is to calculate the cluster form factor since in a percolating lattice 
clusters have different sizes and shapes. We observe however, that if we are sufficiently close 
to the ordering vector q = Q we can write 

F N (q) = e ~R 2 W^-Q) 2 -Q 2 \ (6.29) 

where R{N) is the gyration radius of the cluster and is roughly given by N l l D a (we normal- 
ized the form factor such that F^(q = 0) = 1). Substituting ( |6.29| ) and Q6.12p into (|6.28|) 
and averaging the result with (|6.6|) we easily find 



[ln^o/M] 1 9 tanh((3u/2) , 2/D 2 2 2 -> 

^ {X(q, £*>)} oc (u/u y- x 6XP ^ i ln ( 2 ^o/^)/7j 1 [(q-Q) ~Q \a \ (6-30) 

where c is a constant of order of unit. Notice that for q = we obtain (|6.20|) . Also notice 
that J dqQ {x(q, cu)} ~ S{x(q = 0, u;)} apart from logarithmic corrections. 

A first direct consequence of the Griffiths phase is the weak interplay between q and uj 
which are essentially decoupled. The scattering is peaked around the ordering wave-vector 
Q and the width of the peak is given by 

1 

a [ln(2^ /M/7] 

which is independent of temperature and weakly dependent on the frequency. Actually only 
measurements over many decades of frequency can detect any significant changes in the 
width of the peak. Notice, moreover, that the width of the peak does not depend on the 
distance from the QCP which is given by A oc £ D (p) and therefore we predict the width 
of the peak to be essentially independent on the concentration. All these results, however, 
have a limitation which is the assumption of the gaussian form factor given in ( |6.29|) . 



W « — — ii/D ( 6 - 31 ) 



3. Specific Heat 



An important function in the context of NFL physics is the low temperature specific 
heat which is given by 
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C V (H, T) oc (3 2 dAA x -\E 2 H + A 2 )sech 2 ( (3yE H + A 2 



cx/r A / 

which has the scaling form: 



2 , x 3 sech 2 (x) 
ax- 



(x 2 - (/3^) 2 ) 1_A/2 



In 



(3u 



x 2 - (/3^) 2 



(6.32) 



C v (H,T) = T x g x 



H T 



(6.33) 



Again we will consider only the case of T, Eh -C uj . At very low fields, Eh "C T", we see 
that ( |6.32D reduces to: 



C V (H,T) oc T A (ln/?o;o 



,1-6 



(6.34) 



Observe that in this case the specific heat coefficient diverges at low temperatures: 

C V (H,T) 1 



7C0 



T 



oc ^-^(ln^o) 1 " 



(6.35) 



with the same leading exponent as the susceptibility. Notice that the scaling function is 
such that g\(x — > 0) — > 1. 

On the other hand at high fields, Eh T, we have as leading behavior: 



C V (H,T) oc 



{(3H) 



2+A/2 



(/M)) A 
JJ2+X/2 



-0H 



y2- 



H 



(6.36) 



which has a Sckotty anomaly due to the field. All these results are summarized in Fig.|8|. 
These results are all valid for temperatures above T*. Below T* we expect the crossover 
discussed in the previous section and a deviation from power law behavior. 



VII. CONCLUSIONS 

The problem of the effect of disorder in magnetic phase transitions is probably one of 
the most important problems of modern condensed matter. The intrinsic complexity of 
this problem is related to the fact that at very low temperatures close to a QCP statistical 
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fluctuations due to disorder and quantum fluctuations due to quantum mechanics are in- 
trinsically linked to each other. This is especially true for the alloys of U and Ce because of 
the existence of the Kondo effect. 

We have presented a unified picture of Kondo hole and ligand systems and the role 
played by alloying in the physical properties of these systems. Starting from the Kondo 
lattice model we have shown that it is possible to write down an effective quantum action 
( p.8|) which contains both the RKKY and the Kondo effect by tracing energy degrees of 
freedom outside of a shell of momentum A around the Fermi surface. This treatment allows 
for a direct comparison between the effects of RKKY and Kondo and reproduces, at the 
Hamiltonian level, the famous Doniach argument. We have shown, however, that alloying 
leads to a broad distribution of exchange couplings and therefore to a very inhomogeneous 
environment. This should be contrasted with old approaches to this problem where the 
concept of chemical pressure has been introduced. Just briefly, this concept proposes that 
chemical substitution is equivalent to applying pressure to the system: if one substitutes 
a smaller atom by a larger one then the lattice expands and it is equivalent to negative 
pressure; if we substitute a larger atom by a smaller one the lattice contracts and therefore 
it is equivalent to positive pressure. The next step on the chemical pressure argument is to 
assert that the exchange J over the entire lattice changes uniformly under pressure and the 
system can undergo a phase transition from an ordered state to a Kondo compensated state 
as shown by T/v in Fig.|6|. From our point of view chemical pressure and applied pressure 
have different effects on the physics of these systems. 

A simple example of the large difference between pressure and chemical substitution can 
be captured in specific heat measurements in CeCug under pressure [[73|] and with substi- 



tution of Ce by Th in Cei-^Th^Cug [0. From the experimental data on CeCug one finds 
that the inverse of 7 can be fitted by a linear relation with the applied pressure, P: 

''" 1 1 1.017—^- (7.1) 



l{P) ' Pn 

where 70 = 1.67 J mole -1 K -2 is the zero pressure result and P ma x — 8.8 kbar is the maximum 
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applied pressure in the experiment. Moreover, the unit cell volume in pure CeCu6 at zero 
applied pressure is Vo = 420. 6A 3 while at maximum applied pressure, P max , of 8.8 kbar 
the volume is 416. 5A 3 . Assuming a linear relation between the volume and the pressure (a 
constant bulk modulus) one finds, 

"0 r max 

where (3 = 9.7 x 1CT 3 . In the case of the doped samples it is experimentally established that 
pure ThCu6 has a unit cell volume of 417. lA 3 . Since the contraction of the lattice seems to 
be homogeneous and isotropic with doping one finds that the volume as a function of doping 
can be written as, 

^ = l-« (7.3) 

where a = 8.3 x 10~ 3 . The main idea behind the "chemical pressure" theory is to assert that 
([7. 2|) and ( |7.3|) are equivalent to each other, that is, there is a one-to-one correspondence 
between an applied pressure and doping the system. If one assumes that this is indeed the 
case we can find the correspondence between applied pressure and doping: 

Pa . A . 

-c (7.4) 



1 + Cpc (7.5) 



Pmax $ 

Using ( |7.4|) together with ( |7. 1| ) one finds that the "chemical pressure" theory yields to the 
following relation, 

7o 
.7(c) j P 

where (p = 0.858. Now one has to compare (|7.5|) to the actual dependence of 7 on doping. 
Since we are interested in the low doping regime one uses the experimental data up to c = 0.4 
which can be well fitted by a straight line of the form , 

-7^ = 1 + Cdc (7.6) 

7(c) 

where (d = 0.249. We must stress that ( [7.6| ) is obtained from an extrapolation of the data 
to the zero temperature limit. A more careful look at the data, however, seems to indicate a 
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divergence of the specific heat exponent in Cei- x Th x C\ie. We will however, insist using the 
extrapolation just to check the identity between chemical pressure and real pressure. We see 



the clear difference between the actual doped sample (|7.6|) and the calculation based on the 
"chemical pressure" argument, ( |7.5| ), for the same value of doping. We find that Cp/Cd ~ 3.4 
and therefore a quite large difference between pressure and doping. The conclusion is that 
there is more than a change in the lattice constants or the volume of the unit cell when a 
system is doped. Actually, from our arguments we actually expect that the ration (d / Cp to 
be temperature dependent! 

In this paper we have proposed that the cut-off A which determines the energetic 
crossover between Kondo effect and RKKY coupling is directly related to the number of 
compensated moments in the system which is given by (|3.16|) . Thus, A has to be deter- 
mined in a self-consistent way from the disorder and characteristic energy scales in the 
problem. We argue that this scale is not determined by the Kondo temperature alone but 
also by the RKKY coupling between moments. Since the RKKY coupling involves states 
well deep into the Fermi surface A can be rather large and of the order of the Fermi momen- 
tum. As a secondary consequence of this discussion we show that the exhaustion paradox 
proposed by Nozieres does not apply in the concentrated limit of the Kondo lattice because 
it fails to take into account the high energy degrees of freedom associated with the RKKY 
interaction. 

Within the percolation theory of the Kondo lattice we show that close to the QCP 
strong response from quantum fluctuation magnetic clusters gives singular contributions 
at low temperatures which are known to be Griffiths-McCoy singularities. We have pre- 
sented a phenomenological approach for this problem which allows the calculation of the 



various physical quantities of interest such as magnetic susceptibilities (as given in ( |5.18| ) 
and ( |6 . 1 7| ) ) , dynamical form factors for neutron scattering (given in ( |6.30|) ) and thermo- 
dynamic quantities (such as the specific heat in ( |5.23| ) and ( |6.34j )). The integrals involved 
are rather complicated and sometimes one has to perform them numerically. Moreover, we 
have shown that when a magnetic field is applied to the system there is a recovery of the 
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ordinary behavior due to the quenching of the clusters in the presence of the magnetic field 
(see ( |6.36| )). We have also discussed the various scalings in the presence of a magnetic field 
which are given in ( |6.15[ ) and ( |6.33| ). 



We have also reviewed the single impurity Kondo problem and have shown that in the 
anisotropic case, which is relevant for the case of spin-orbit coupling, it can be interpreted 
as the dissipative dynamics of a two level system in a transverse field. We have shown that 
while the XY component of the Kondo coupling acts as a transverse field on the magnetic 
spin, the Ising component gives rise to dissipation due to the production of particle-hole 
excitations at the Fermi surface. Using exact results for the dissipative two level system 
we have shown that the physical quantities differ slightly from the non-dissipative two level 
system by the presence of the parameter a that controls the coupling of the two level system 
to the electronic heat bath. We have also revisited the two impurity Kondo problem and have 
shown that the stable fixed points are related to the coherent flipping motion of the spins 
as a single degree of freedom (as shown in ( f4.54Q and ( |4.58| )). This is a clear sign of cluster 



behavior, that is, coherent quantum mechanical motion. In a cluster of N atoms we have 
identified two mechanisms for quantum tunneling: in XYZ magnets the anisotropy of the 
RKKY interaction generated coherent cluster tunneling (with a cluster Kondo temperature 
given in ( [4.77D ) while in XXZ or Heisenberg systems the Kondo effect generates the tunneling 



(with a cluster Kondo temperature given in ( |4.86| ) ) . In both cases the tunneling splitting 



is an exponential function of the number of atoms in the cluster as shown in ([4 . T0|) . This 
calculation allows us to relate the microscopic quantities in the Hamiltonian (|3.8|) and the 
phenomenological parameters in Section |V]]. We have shown that there is a cluster Kondo 
effect which is parameterized also by a dissipative constant a(N) = (N/N^ where ip is 
an exponent which depends on the type of coupling with the particle-hole continuum (see 
(TO)) and N c (given in (^83|)-(^85|) and (^88|)-(^89|)) is the largest size of the cluster 



for which the Kondo effect ceases to exist and the cluster are frozen. Thus, the Kondo 
temperature of the cluster vanishes at N c . In this respect our picture is quite close to the 



Kondo disorder picture proposed for the explanation of NFL in these systems |Tl| , |21|| . We 
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have to point out, however, that the microscopic origin of both pictures is quite different 
since in our scenario the system has to be close to the magnetic phase transition and the 
Kondo temperature is defined for an extensive object which is determined by the laws of 
percolation theory. 

Based on these results we proposed a dissipative droplet model for the magnetic behavior 



in disordered metallic magnets which is described by the action (|5.4f ). This picture is the 



extension of the quantum droplet model proposed by Thill and Huse Jf2| for the case of 
insulating magnets. We have shown that in this model there is a crossover temperature 
T* « Epe~' yNc '' fi above which the results of the Griffiths-McCoy phase are still valid with 
power law behavior in many physical quantities (see (|5.12|) ). Below T* the behavior of 
the physical quantities is modified by dissipation and a new behavior emerges where the 
magnetic susceptibility and specific heat diverge stronger than power law and do not scale 
together as shown in ( |5.18|) and ( |5.23|) . Actually the magnetic susceptibility has a stronger 
divergence than the specific heat coefficient. This could perhaps explain why the exponents 
A found from magnetic susceptibility measurements are systematically smaller than the ones 
found in the specific heat data pp|J75[ ]. We have estimated, however, that T* is rather small. 
Thus, power laws as given in ( |6.17| ) and fl6.34| ) should dominate the crossover behavior at 



low temperatures. We have calculated the thermodynamic properties of the clusters alone 
but it is clear that in a real system one has to add also the contributions coming from the 
paramagnetic environment (Fermi liquid or heavy Fermi liquid) which are given by ( [2.1|) and 
( p.4| ), for instance. Since the response functions we have calculated are singular as T — > 
they dominate the response at low enough temperatures. 

Finally, we should mention that our picture is based on two main features: the short 
range character of the disorder and the complete decoupling of the clusters. On the one hand, 
if the disorder is long range correlated then stronger divergences of the physical quantities 
are possible |76j . On the other hand, if the cluster are not decoupled but there is a residual 
interaction between them it is possible that the system reaches a spin-glass state at lower 
temperatures. The nature of the spin glass state depends on the dissipation. If the freezing 
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temperature is of order of T* then the dissipation will kill quantum coherence and a classical 
spin glass state is formed. If T* is smaller than the freezing temperature then a quantum 
spin glass ground state is possible |72| . The existence of such phases can only be decided on 
experimental basis. 

Thus, in this paper we put forward the theoretical basis for Griffiths-McCoy singularities 
in the paramagnetic region close to the QCP for magnetic order of U and Ce intermetallics. 
We show that strong temperature dependence in the physical quantities is expected at low 
temperatures due to the quantum mechanical response of magnetic clusters. We believe that 



our picture is in agreement with a large number of experiments in these systems [[j^.[39|.|75|] . It 
would be interesting to investigate further the existence of a even lower temperature regime 
where dissipation plays a decisive role in the physics of these systems. 
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APPENDIX A: PERTURBATION THEORY 



Our starting point is (|3.7|) where we integrate out the electrons in regions f^i and and 
obtain an effective action for the electrons residing in region Q : 

Z = J DS(n,t)D%{r,t)DMr,t)e iSeff(s ^ oM (Al) 

where 

S ef , = S - itr ln(G _1 ) - z £ j dt J dt ! £ 5X t (r, t)G&(r, t; r', *>av(r', t') (A2) 

r,r' ' a,a' 1,1! 

where Sq is ( |3.4j) for the electrons in region Q only. Furthermore, 

riaAr, *) = £ J a , b (r)S a (r, t)r b a ^ {r, t) , (A3) 

a,b 

and 

(r, t, l, alG' 1 ^', t', if, a) = + fi + e a (-iV)J 8 rir >8(t - t')5^5 a ^ 

+<' E Ja, b (r)S a (r, t)rl p 5{t - t')8 r y (A4) 

a, 6 

where erf^ is a Pauli matrix which works in the subspace of states Oi and S7 2 - Moreover, 

4,fc(r) = E ^a,b(«)5 r ,r n • (A5) 
n 

It is convenient to define the bare Green's function as the solution of the J = problem, 
that is, 

(i^ + li + eai-iV^j G 0yCt (r -r',t- t') = 8(t - t')8 IjT , . (A6) 

Notice that the time ordered Green's function can be trivially obtained from the above 
equation and reads, 

Go,a(k, u) = — — ■ r (A7) 

to + fi - e a (k) - -sgn(e a - fj,) 

where we have introduced the electron relaxation time r due to weak disorder. Thus, in the 
regions of momentum space introduced in Section ( p|) we have 
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Ga , M - e <*'-- A -*> 



8(fc - k F , a - A) 

G Q)a , 2 (k,a;)- a; + ^_ ea(k) _ i . (A8) 



Notice that the trace over the fast modes has generated two kinds of terms in (IA4): new 
interactions between the flo electrons and the spins (through if) and spin-spin interactions 
through G. In what follows we will consider perturbation theory of J for these fast modes. 
First let us consider the new interaction terms between electrons and spins. Since rj is 
already order J we can consider G ~ G in ( [A"^ ) and rewrite 



I = J2[dt[ dt'EEW r >*) G X>>*; r ^W(*V) 



a, a' 



/ = E/^/^'E^,oM) 

r,r' 7,5 



E EE^W^(r')^(r,t)5 c (r',t') 

a,b,c,d a,/3 

r^G^ a! (r,t-,r',t')}^,o(r',t'). (A9) 



Thus, in order to evaluate / one needs an expression for G. From ( |A4[ ) we have (symboli- 
cally), 

G- 1 = G 1 + J = G 1 (1 + G J) 
G = (l + G jy 1 G 

« G + G 0l 7G + /i.oi. (A10) 

where J7 is the second term in ( |A4| ) . 
To leading order one has, 

I^Y.j dt j dt'Y^^M^'J) E JaAr)JcAr')S a (r,t)S c (r',t') 

r.r' 7,5 a,b,c,d 

x E<7<* (E G <w(r-r',t-t')) (All) 



which is a electron spin-flip scattering process involving two different local moments. Observe 
that this interaction is mediated by the single particle Green's function Go which can be 
evaluated from ( |A8l) . After a trivial frequency integral we find 

77 



np-\ t \l T r rile 

E<Wr,t) = J (^s^-(*-^) [9(fc - - A)0(t) 

-e(^, a -A-fc)e(-t)] (A12) 

where pi = N/V is the lattice density. Since we are only interested in the asymptotic 
behavior of the system we observe that when t — > oo the integral is dominated by the saddle 
point at /i = e a (kF,a) but the functions are only finite outside this region, thus, the above 
integral has very fast oscillations given null contribution. In order to see this result explicitly 
let us expand around the saddle point, that is, we change variables: 

k = kp, a + Q u F,a (A13) 

where uf jQ is a vector perpendicular to the Fermi surface (that is, in the direction of the 
Fermi velocity). Thus we can write 

e k « n + v F>a q 

k^k F + q (A14) 

and the integral above becomes 

E GW(r, t) « E TTT^ / dge^'--^'.-)-'-^'-*] [6(g - A)0(t) 

PL kFQ (2tt) 3 J-oc 

-e(-A-?)e(-t)]e-* (A15) 

where we introduced an infinitesimal converging factor v — > 0. The integrals are now trivial 
and give 

l2 ~-*A(u J ? l „-r-t> Fo[ t)aflw(t) 

E GW(r, t) « E ^^" r - ~ ~ 7 -—77^ ( A16 ) 



k 



F.a 



7r 3 pi uj7 ]Q • r — t>ir at — iusgn{t) 



which shows that at long times the Green's function oscillates with frequency v F A and 
therefore averages out to zero. 

We now look at the the first term in (|A2|) and write 

tr ln(G 1 + J) = tr ln[G \l + G J)} = trlnGg 1 + *rln(l + G J) 

« tr In Go 1 + tr(G J) - ^tr(G JG J) + h.o.t. (A17) 



In order to simplify the notation we will use a simple dummy index % = (r, t, a, t) for the 
matrix elements. Observe that tr{GoJ) = ^W^j <Ji,j — 12iGo,i^Ji,i = because of = 
in (|A4). The second order term in ( |A17| ) becomes 



tr(GoJ GqJ) — X Go,i8i,jiJj,kGo,k8k,i<?k,i 

i,j,k,l 

= X/ X/ Go,i>?i,kGo,k»?k,i • (A18) 
i,k a,7 

Using ( |A4|) we find explicitly 

tr(G JG J) = Jdrjdr' JdtJ dt' X -Mr) J c , d (r')S a (r, t)S c (r', t')U b4 (r - r', t - t') (A19) 

a.b.c.d 

where 

n M (r,t) =2EE<7i G o,«. i (M)< 1 -Go, 7 , ( fr, -f) (A20) 
is the polarization function of the electron band. It can be rewritten in momentum space as 
n M (p,^) =JEEW/^ / ^ / 7^ G o w (k^)Go, 7/ (k + p^ + fi) (A21) 

which can be evaluated once the conduction band is known. Thus, to second order the 
integration of the fast degree of freedom generates an interaction among the spins in different 
sites. One further approximation which is usually used is to take the static limit in (|A19|). 



This approximation is good if one looks at times scales which are much larger than the 
scattering time of the electrons between moments which is £/vp ja where I is the distance 
between moments and vp >a is the Fermi velocity for the a branch. We replace n& it j(p, u) 
by K[n bd (p, 0)] (since S[II bjd (p, 0)] = 0). In this limit the second order contribution is the 
RKKY interaction between localized spins (which depends on the cut-off A) which can be 
written as 

H mKY (A) = X X T a,b(rn - r m , A)S a (r n )S b (r m ) (A22) 

m,n a ,b 

where 

r a>6 (r, A) = - £ J a , c J M j rfqe iq - r 3?[n c , d (q, Q = 0, A)] . (A23) 



c.d 
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The integral in ( |A23| ) can be written directly from (|A21 ) and reads 



dk! I 7- ; , e 



+ 



i(k-k')-r^ / Q(fcF,« - A - fc)9(fc' - k Fa - A) 
(27r) d ~ 2tt \(lu + - e a (k) - i/r)(u + pi - e 7 (k') + i/r) 

0(fc - fc F ,« - A)6(A; Fi7 - A - k') \ 



{u + pi — e a (k) + i/r){uj + pi — Sy(k') - i/r) 
dk 



i / dk' 



i(k-k')-r [ ®( k F,a - A - k)&(k' - k Fn - A) 



+ 



(2vr) rf \ e 7 (fc') - e a (Jfe) + 2i/r 

Q(k - k F , a - A)6(fc F , 7 — A — k')\ 



(A24) 



e a (fc) -e 7 (fc') + 2z/r J 

Let us consider the special case where the dispersion are the same in both spin branches 
and the Fermi surface is assumed to be spherical. In this case the bare Green's function 
does not depend on the spin label and n 6j(i (r, t) = 5f,,dII (r, i). We can write (|A24 ) as 



J(r, A) = i I dk' 
+ 



dk 



3 i(k-k') 



k ,v r (Q(k F - A - k)Q(k' - k F - A) 



(2tt) 3 \ e(Jfe') - e{k) + 2%/t 

3{k — k F — A)Q{k F - A - k') \ 



dk' 



e(k) - e(k') + 2i/r / 
cos[(k - k') • r] 
(2tt)3 e(Fy^(A?)T2i7r 



Q(k F - A - fc^A/ - k F - A) 



2 rk F -A roc kk'sm(kr) sin(fc'r) 
nr 2 Jo Jk F +A e(k') — e(&) + 2i/r 



(A25) 



When A — > and ^> 1/r we find the usual result [52 



Z(r,0) = ^e^F (2A*r) 



(A26) 



where 



k F r 



m 



v F r 



(A27) 



is the electron mean-free path and 



7T 



TAx) = - 



sin(x) — x cos(x) 



x/- 



(A28) 



which leads to the final expression for the RKKY interaction in d = 3. 
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Let us now consider the corrections to the RKKY interaction due to a finite A/fcp. In 
all the cases here we will consider the situation in which kp > A ^> 1/r. From ( |A25| ) we see 
that the correction has the form 



4m k 



51(r, A) = 7 fe- r/ "J r (2k F r, A/k F ) (A29) 

7ir z 



where, T = T\ + JF 2 + JF 3 , with 

~ x , r 1+ v , zw sinixz) sin(xw) 



F x (2x,y) = - J dz J dw _ 2 



T%({Lx,y) — — dz dw 

Jl—V Jl 



Z z — W z + IT] 

00 , zw sm(xz) sin(xw) 



i-y Ji z 2 — w 2 + irj 



F 3 {2x,y)= I dz I dw 2 ^ ; ; t J (A30) 



ri+y zw sin(xz) sm(xw) 
1 dw 2 ^T - ' 

'l-j/ Jl Z A — W z + 17] 

where 77 — ^ 0. All the integrals can be evaluated after a very tedious algebra. The final 
result reads 

!F{x,y) = -|— ^ |vrcos(x) — (n — 2xy) sin(x) — 2xsin(xy) — 2x 2 y{Ci{x) — Ci(xy)) 
+ cos[(l + y)x] [(Ci(x) — Ci(xy)) — (1 + y)x(Si(x) — Si(xy))] 
+ sin[(l + y)x] [(Si(x) — Si(xy)) + (1 + y)x(Ci(x) — Ci(xy))} 

— cos[(l — y)x\ [(Ci(x) — Ci(xy)) — (1 — y)x(Si(x) + Si(xy))] 

- sin[(l - y)x\ [(Si(x) + Si(xy)) + (1 - y)x(Ci(x) - Ct(xy))}} (A31) 

where 

Ci(x) = - f dt^ 
» rft sin(t) 

X 



Si( x ) = - [°° dt 8 ^- (A32) 

J X t 



are cosine and sine integrals. A plot of ^(x, y) is shown in Fig.[| Thus, the actual RKKY 
interaction in the presence of a finite A/ kp is given by the sum of ( |A26| ) plus ( |A29| ) . Observe 



that besides the usual 2kp oscillations a finite A/kp produces other oscillations at 2A. Thus, 
the structure of the RKKY interaction is more complicated. 
At very short distances, 2kpr <C 1, we have 



X(r,A) 



2mkp 



3r 




(A33) 



which, as in the case of the usual RKKY interaction, leads to ferromagnetic coupling. The 
main difference, however is that the slope of the RKKY is reduced by the factor (1 — 
A/ hp) 3 and therefore the strength of the ferromagnetic coupling is reduced as A increases. 
Furthermore, while the first zero of the RKKY interaction for A = occurs for kpr pa 2.2467; 
when A ^ the first zero is reduced as shown in FigJ4| and saturates at kpr pa 1 when 
A pa hp. Overall the effect of a finite A is to increase the range of the antiferromagnetic 
coupling for small to moderate kpr. The component near the Fermi energy contributes a 
largely ferromagnetic component at these kpr. Note that overall the contributions to RKKY 
from the energy region away from the Fermi surface is at least comparable to that for near 
the Fermi surface. At intermediate distances, that is, A -1 ^> r > kp 1 , the RKKY interaction 
can be written as 



Til nj TP i\. 

J(r, A) pa l(r, 0) H ^— sin(fc F r) [-tt cos(k F r) - 2(1 - C + Ci(2k F r) - ln(2Ar)) sin(A; F r) 



where C ~ 0.577216 is the Euler constant. We see from these expression that the correction 
to the A = case vanishes as (A/kp) hx(A/kp) as A/kp — > 0. The most striking difference 
between the behavior of the RKKY interaction with A/kp ^ is it behavior at large 
distances. It is very simple to show that for r » kp 1 , A -1 one has 



which decays like 1/r 4 instead of the usual 1/r 3 . Thus, for finite A the RKKY interaction 
has a shorter range. 

APPENDIX B: MEAN FIELD THEORY OF ITINERANT FERROMAGNETISM 

IN THE KONDO LATTICE 



+ 2cos(k F r){Si{2k F r)+7T/2)) 



(A34) 




(A35) 



In this appendix we will discuss the physics of Hmf in ( |3.13| ). Assuming that the 



magnetic order is homogeneous we can write 
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N s N a 



Hmf = E e<x( k )Ck, CT Ck, CT + J Z M Z ^(n iiT - n id ) + J z m z ^ S* (Bl) 

k.cr i=l i=l 



where 



"•^jE^-W) (B2) 

where iV, is the number of magnetic atoms. As usual we can now perform a Fourier transform 
of the electron operator 

V = 4E^" c M (B3) 
V Jv k 

where N is the total number of atoms. We rewrite (IB II) as 



Hmf = E [( e t( k ) + J z cM z ) c k)T c k , T + (e ; (k) - J z cM z ) c ka c k , 

k 

+ J z m z J2S* (B4) 

j 

which brings the Hamiltonian to diagonal form (c = N s /N). In the ground state we have 

(n k)CT ) = e(A*, ff - fc) (B5) 



and therefore 



1 ^ 1 

M = TT EK^ = TtE^) 
iV s i=l iV k 



1 fc 3 



(B6) 



PL 67T 2 

where Pl = A^/V is the lattice density. On the other hand, the total number of electrons is 

N 

N e = T, E(V> = E EKa) 
i=l o" k o- 

= ^V% (B7) 
and thus the number of electrons per lattice site is 
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(B8) 



moreover 



m z — n-j- — rij_ (B9) 

from which we conclude 

n + m z 
n-t = 

_ 2 

n — m z 

ni = —^-. (BIO) 

Since the up and down spins are in thermal equilibrium they must have the same chemical 
potential, that is, 

^t-. + J ' M ' c = t>- J ' M ' c (B11 » 



which can be rewritten with the help of ( |B6|) as 



((-T)) 2/3 -((ni)) 2/3 = -^^. (B12) 



Using ( |B1Q ) we find 



(n + m z y/* - (ft - m 2 ) 2 / 3 = (B13) 

which has to be solved for m z . If one is interested in the critical temperature T c of the 

problem which is much smaller than fi it is sufficient to consider the case of m z <C n in 
which case we find 

• <B14) 

On the other hand, from ( (BT| ) we have 

M 2 = -SB S ( P SJ ><™A . (B15) 



The solution of the problem is given by the set ( |B14p and ( |B15p . In particular, substituting 
(|BT1 into (|BT5| ) 
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Mz = SBs ( ^f&T ".) (B16) 



2(3tt 2 Pl ) 2 /3 

At T — > T c ~ we have M 2 — > and therefore we find 

S(S + l)m*J 2 c 2 (n) 1 / 3 

2(37r 2 p L ) 2 / 3 
S(S + l)c 2 nJ 2 S(S+1] 



N(0)J 2 (B17) 



4 £f 6p L 
which agrees with (|3.9|). 

APPENDIX C: FREQUENCY DEPENDENT SUSCEPTIBILITY OF A TWO 

LEVEL SYSTEM 



Let us consider the response of a system described by the Hamiltonian ( |6.7D to an oscil- 
lating magnetic field of frequency u. In linear response this is given by the retarded spin-spin 
correlation function averaged over the thermal ensemble 

X (t) = -iQ(t)([Sa z (0),5a z (t)}) = ~^ tr (e- pH [Sa z ,e iHt 6a z e- iHt )) 

= (e-^» - e~ pEm ) e^-^Kn^lm)] 2 (CI) 



run 



where H\n) = E n \n) and Z = J2n e n an d = &z — (&z)- A simple Fourier transform 
gives 

\{n\5a z \m)\ 2 



-f3E„ 



z V J uj + (E n - E m ) - it 

4s( 



e -p En _ e -p Em \ \(n\6a z \m)\ 2 ( _ 



Z n,m E n — E m \ lu + (E n — E m ) —iej 

oj / _p E _ j3E \ I (n\5a z \m) \ 2 , 

= X[ ) + ^h [e m) (E n -E m )(u + (E n -E m )-ie) (C2) 

where e — f and 

« = 4E(e*-«*)*f (03) 

n,m 

is the static susceptibility. 
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The calculation is very simple because we are dealing with a two-level system problem. 
Indeed the Hamiltonian ( |6.7| ) can be diagonalized by a unitary transformation 



U = e l0tT v 



tan(20) 



E 



(C4) 



H 



which transforms 



UHcU- 1 = J El + A 2 t, 



UaM 



-i 



1 



Ua x U 



-i 



Ejj + A 2 



E 2 H + A 2 



[E h t z - Ar, 



[At, + £ h t x ] . 



(C5) 



Moreover, in the new basis of the Pauli matrices r with eigenstates |±) and eigenenergies 



±J E H + A 2 , respectively, we can easily show that 



(±KI±> 
(±k*lT> 



tanh ^/^ + A 2 
^ + A 2 V^ V 



^ + A 2 

_A 

E 2 H + A 2 



(C6) 



Substituting ( |C6|) into ( P3|) we easily find 

(3E H 



x(o) 



E?r + A 2 ' 



sec/i 2 + A 2 + 



A 2 



tanh + A 



(C7) 



^T" , , (^ + A 2 )3/ 2 

which is the result ( ft. lip . The dynamic part can be obtained in analogous way and it gives 



xM = x(o) 



E]j + A 2 ' 



sech' PJE 2 H + A 



a; 



w — it 



(^ + l 2 )3/ 2 tanh W^ + A2 ) 



+ 



_cu - 2y^| + A 2 - it u + 2y^| + A 2 - ie 
from which we can extract the real and imaginary parts: 
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a [xH] = x(o) 



E 2 H + A< 
A 2 



sec/i^ [pJEl + A* 



oo 2 + e 2 



2{E 2 H + A 2 ) 3 / 2 



tanh pJE 2 H + A 



w(w - 2JE% + A 2 ) 



- 2JE 2 H + A 2 ) 2 + e 2 



+ 



tufa; + 2 J Elf + A 2 ) 



u + 2^E 2 H + A 2 ) 2 + e 2 
A 2 



+ A 2 

which are used in the calculations. 



tanh /3JE% + A 2 



6(u - 2JE 2 H + A 2 ) - 5(u + 2JE% + A 2 ) 



(C9) 



APPENDIX D: TWO IMPURITY CALCULATIONS 

Starting from ( |4.36| ) and using the expansion ( |4.37| )the Hamiltonian is written as 
dk 



2 



f dk 


r dk' r 


1 2^K J 


2n - 



a ,j,l,a,/3 



- JmAk, k')^k)rlK^iAk'){S l>a - S 2 , a ) + J-, a (k, ^Jk^a^M^S^ + S 2 , a ) 

- iJ t>a (k, k')^lAk)rlK,^iAk')(Si, a - S 2> a)} (Dl) 

where r a are Pauli matrices which act on the subspace of j = 1, 2, and 

kk' T R 

J±,a = (N e (k)N e (k') ± N (k)N (k')) 

kk' T R 

Jm(t), a = — -2- (N e (k)N (k r ) ± N (k)N e (k')) (D2) 

are the exchange constants. Observe that all the momenta here are defined in a thin shell 
around the Fermi surface. Thus, we can linearize the band by writing = vp{kp — k) and 
expand the above exchange constants to leading order in k in which case we get 



R 



J +ta (k,k')aicN(p)J2 

j.,( fc>ifc 0«^(o)j?^^ 
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J m , a (k,k')^TTN(0)jK 

J t , a (k,k')^0. 



\ 



sin(/ci?i?) 
k F R 



(D3) 



in which case there is no momentum dependence in the coupling constants of ( pi| ) which 
transforms it to an impurity problem if we define the Fourier transform: 

"+°° dk 



/+oo 
-oo ZTT 



(D4) 



in which case ( |D1| ) becomes turns into (|4.39|) . 

An important term we have not included is due to ordinary scattering of electrons at the 
impurity. This term is just the electrostatic coupling which is given by 

H v = ]T (y(R/2)V+ (R/2)^«(R/2) + V(-R/2)V£(-R/2)V> a (-R/2; 

a 

= E E (v(K/2)e-^-^ 2 + V(-R/2)e^ k - k ')- R / 2 ) c^q,, Q 

a k,k' 



(D5) 



/• dk 


r dk' r 


1 2~H ) 


2^T - 



where we used ( [4.37] ) . As previously, assuming that the electron momenta is close to the 
Fermi surface we can rewrite the above term as 

H v = E [Vi4,a( x = 0)^M* = 0) + V 2 4, a (x = 0)il> 1>a (x = 0) 

a 

+ V (4A X = tyfoA* = 0) + 4jx = 0)ipx, a (x = 0))] . (D6) 

The bosonization procedure goes like the one described for the single impurity problem 
in Subsection [IV A| . The only difference is that we are going to construct the factors K a 
explicitly in terms of the fermion fields. We define 



i^2,cr = exp linjdx 



(D7) 



which can be rewritten in terms of the bosonic fields if we use that 



1 d^ a 



(D8) 



2n dx 

It is easy to show that ( |4.39|) reads (for simplicity we work again with the problem of uniaxial 



symmetry) 



H 



+ 



^ i=c,s,f,sf 

v F J +jZ <9$ s (0) 



dx 



2V 



+ r a Sl,aS 2 ,a + — COS($ S/ (0)) COS($/(0) 



na 



2n dx 

+ !^±^ cos ($ s/ (o)) (cos($ s (0))(S li:E + S 2>x ) - fflB($ a (0))(5 liW + 5 2iW )) 



+ 



7ra 

27r dx 



r/ / " 1 sin($ s/ (0)) (sin($ s (0))(S 1 ,, - 5 V ) + cos($ s {0))(S 1>y - S^)) 



7ra 

+ V -I±L sin($ s/ (0)) sin^O) - 0)(S 1>t + S 2> 



+ 



v fJ-, 



cos($ / (0) - 9) (cos($ s (0))(S M + S 2>x ) - sin($ s (0))(5 liS/ + S 2>v )) 



(D9) 



which can be rewritten in a more economical format: 

/ o , \ 2" 



2 . 



i=c,s,f,sf 



U(x) 



dx 



+ 2J F a Si ta S 2ta 



+ — cos($ s/ (O))cos($ / (O)-0) 



7ra 



v F J+, z d$ s (0) , c 
+ ^ &-^ + ^ 

+ cos($ s/ (0)) (e"-W(5i, + + S 2 , + ) + e-w-^CV + S 2 ,_)) 



+ 



27ra 

w F J m , 2 d$ s/ (0) 



27r cte 



($ s/ (0))(e^( )(5 li+ + S 2 ,_ 



— e 



+ 2 sin 

2na 

+ ^Zii ^(^(0) - 0) sm($ s/ (0))(£ M + S 2>z 



-i*«,(0) 



(St,- + s 2 , 



+ 



na 
VfJ-,. 
2na 



cos 



($/(0) - 0) (e^ (0) (5 li+ + 5 2 ,+) + e- i *'^(S' 1 ,_ + &,-)) 



(D10) 



where S it ± = S i>x ± iS^ y and 6* is defined in ( |4.44| ). Notice that the charge degrees of 
freedom have decoupled entirely and can be disregarded and that the diagonal terms in 
can be written as gradients of the boson operators and therefore can be trivially 
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absorbed into a shift of the bosons and do not appear here. After the unitary transformation: 



U 



s +Sa,*)*.(o) the Hamiltonian becomes ( |4.43|) . 
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FIGURES 

T(x) /T(0) 

c c 

Paramagnetic 
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FIG. 1. Typical phase diagram of the systems studied with an ordered magnetic phase and an 
anomalous metallic phase close to the QCP. The Griffiths-McCoy singularities appear close to the 
QCP. 
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FIG. 2. Fermi Surface and its regions of ener; 
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FIG. 3. Plot of J-(2x,y) which appears in the expression for the RKKY interaction ( |3.9|) : 
dashed line is the usual RKKY (A = 0); continuous line A/Jzf = 0.1. 
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FIG. 4. Plot of the first zero of the RKKY interaction, y = kpro, as a function of x = K/Uf 
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(b) 




FIG. 5. Percolation problem in a (a) Kondo hole system and (b) Ligand system. p e is the elec- 
tronic density; pf is the density of magnetic moments; pq is the density of compensated moments; 
p is the percolation parameter and p c is the percolation threshold. 



Energy Scales 



RKKY 




RKKY Kondo 

FIG. 6. Doniach phase diagram: long dashed line is the Kondo temperature, Tjc; short dashed 



line is the RKKY temperature,TR^-^y; the continuous line is the ordering temperature TV. 
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\w 
— 

iTtF 
\w 



12, -> 

= <* * J 2 /f 

12, +> 



r/2 

z 



II, +> 

= < ; j 2 /r 
li, -> 

FIG. 7. Energy level splitting in the two- impurity Kondo effect. 



.X 



g^H/T) ~1 
f (H/T)~l 
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\ 

/* f ^( H/ T ) ~ ( T /H ) ^ \ 

2+X/2 -H/T \ 

/ g (H/T) ~ (H/T) e 1 



H 
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FIG. 8. Behavior of the scaling functions of the magnetization, f\(H/T), given in flOp and 
the specific heat, g\(H/T), given by ( |6.33 ), in the H x T phase diagram. The line T ~ H marks 
the crossover in behavior and T is the value of the RKKY interaction above which the cluster does 
not exist. 
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